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Harmonic  Gyrotron  Oscillators  and  Amplifiers 


1  Summary 

This  report  summarizes  the  work  done  on  the  task  titled  “Harmonic  Gyrotron  Oscillators  and 
Amplifiers”,  covering  the  period  from  September  30,  1994  to  September  29,  1995.  During  the 
past  year  we  have  completed  the  test  experiment  of  phase  locking  a  35  GHz  second  harmonic 
gyrotron  oscillator  using  a  quasi-optical  circulator  to  separate  injection  and  output  signals. 
The  experiment  well  demonstrated  the  capabilities  of  the  hot  test  equipment  in  the  harmonic 
gyrotron  laboratory  at  UMCP  and  provided  a  good  basis  for  understanding  the  importance  of 
mode  selective  design  and  the  phase-stable  operating  parameters  of  more  advanced  harmonic 
gyrotrons. 

We  have  also  developed  a  special  approach  to  realize  a  compact  high  performance  MMW 
source.  This  special  approach  includes  two  aspects.  First,  we  employ  a  new  gyrotron  op¬ 
eration  regime  in  which  multi-stage  interaction  is  used  with  resonant  cyclotron  harmonic 
multiplication  in  each  succeeding  stage  to  provide  sufficient  prebunching  of  the  MIG  pro¬ 
duced  electron  beam  at  lower  harmonics,  this  will  provide  an  efficient  excitation  of  the  EM 
field  in  the  final  stage  of  the  device  at  a  higher  harmonic  of  the  electron  cyclotron  fre¬ 
quency.  Secondly  we  employ  an  experimentally  demonstrated  TEo„  mode  selective  circuit 
for  input  coupler  and  for  each  interaction  stage  in  order  to  suppress  unwanted  competing 
modes  which  are  inevitable  in  any  harmonic  gyro-device.  For  the  purpose  of  demonstration 
of  this  approach,  two  proof-of-principle  devices  have  been  designed  and  contracted  for  fabri¬ 
cation.  These  devices  are  the  two-stage  harmonic  multiplying  gyro-TWT  and  the  two-stage 
phase- controlled  harmonic  multiplying  inverted  gyrotwystron  (phigtron).  Both  of  them  are 
with  Ku-band,  fundamental,  gyro-TWT  prebunching  and  Ka-band  second  harmonic  output. 
In  addition,  we  have  also  made  good  progress  in  preparing  a  Ka-band  tunable  magnetron 
source  that  might  be  used  for  driving  a  gyro-TWT  with  fourth  harmonic  W-band  output 
and  second  harmonic  Ka-band  prebunching.  The  prebunched  beam  might  also  be  used  to 
optimize  the  efficiency  of  a  fourth  harmonic  oscillator.  This  study  on  harmonic  gyrotron 
amplifiers  and  oscillators  has  also  been  supported  by  MURI  program  since  May  1,  1995.  As 
a  result,  the  research  works  are  being  able  to  continued. 


2  Accomplishments  to  Date 

During  this  reporting  period  and  the  few  months  preceding  it,  significant  progress  was  made 
both  in  the  experimental  and  theoretical  aspects  of  the  program.  Accomplishments  to  date 
include  the  following: 

1.  Completion  of  the  test  experiment  of  phase  locking  a  35  GHz  second  harmonic  gyrotron 
oscillator  using  a  quasi-optical  circulator  to  separate  injection  and  output  signals.  De¬ 
tails  of  the  test  and  analysis  of  the  experimental  I'esults  are  presented  in  the  appendix 
of  this  report. 


2.  Development  of  a  new  gyrotron  concept,  “harmonic  multiplying  gyrotron  amplifiers 
and  oscillators”,  with  the  potential  for  weakly  relativistic  small  orbit  gyro-devices  to 
efficiently  operate  at  fourth  harmonic  (Figures  1  and  2). 

3.  Completion  of  the  engineering  design  of  the  35  GHz,  second  harmonic  gyro-TWT 
amplifier  and  the  inverted  gyrotwystron  oscillator  (phigtron)  both  implemented  with 
complex  mode  selective  circuitry  and  Ku  band  fundamental  prebunching  (Figures  3 
and  4  and  Table  1). 

4.  Contracting  for  fabrication  of  the  above  35  GHz  proof-of-principle  harmonic  multiply¬ 
ing  gyrotrons. 

5.  Procurement  of  a  surplus  radar  pulsar,  16  kV/16  A,  with  reconditioning  and  some 
modification  as  the  high  voltage  modulator  for  a  Ka-band  tunable  magnetron  which 
has  been  borrowed  from  Prof.  K.R.  Chu  (National  Tsing  Hua  University,  Taiwan, 
R.O.C.)  and  is  planned  for  driving  our  fourth  harmonic  gyrotrons.  This  modulator  has 
been  installed,  adjusted  and  made  to  work  adequately  in  the  magnetron  performance 
test  experiment. 


3  Recent  Papers  and  Presentations  Related  to  this 
Task 

1.  H.  Guo,  V.L.Granatstein,  B.  Levush,  T.M.  Antonsen,  Jr.,  J.P.  Tate,  J.  Rodgers,  and 
W.W.  Destler,  “Experimental  Study  of  Phase- Controlled  Harmonic  Gyrotron  Oscil¬ 
lators  and  Amplifiers,”  1994  Microwave  Power  Tube  Conference,  Monterey,  CA,  May 
10-12,  1994,  abstracts  in  conference  record,  presentation  2B.7. 

2.  J.P.  Tate,  H.  Guo,  M.  Naiman,  L.  Chen  and  V.L.  Granatstein,  “Experimental  Proof- 
of-Principle  Results  on  a  Mode-Selective  Input  Coupler  for  Gyrotron  Applications,” 
IEEE  Trans.  MTT  42,  1910  (1994). 

3.  H.  Guo,  V.L.  Granatstein,  T.M.  Antonsen,  Jr.,  B.  Levush,  J.P.  Tate,  and  S.H.  Chen,“ 
Compact,  Harmonic  Multiplying  Gyrotron  Amplifiers,”  1995  IEEE  Int.  Conf.  on 
Plasma  Science,  June  5-8,  1995;  abstract  in  Conference  Record  95-CH35796,  p.  211. 

4.  H.  Guo,  D.J.  Hoppe,  J.  Rodgers,  R.M.  Perez,  J.P.  Tate,  B.L.  Conroy,  V.L.  Granatstein, 
A.  Bhanji,  P.E.  Latham,  G.S.  Nusinovich,  M.  Naiman,  and  S.H.  Chen,  “Phase-Locking 
of  a  Second  Harmonic  Gyrotron  Oscillator  Using  a  Quasi-Optical  Circulator  to  Separate 
Injection  and  Output  Signals,”  IEEE  Trans.  Plasma  Sci.,  Oct.  1995,  to  be  published. 
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Figure  1  Physical  Mechanism  and  Configuration  of  the  Harmonic 
-  Multiplying  Gyrotron 


1st.  stage:  i)  beam  modulation  by  Vf,  ~  J  Ei'^'“‘d(re0) 

.  ii)  fundamental  ECRM  interaction 
iii)  %  =  Iq  entrance 

2nd.  stage:  i)  =  Iq+Ij-j  +  Ij2  ••••  at  the  entrance  and  f2  =  2  fj 

ii)  Vf2  ~  J  E2  d(rj6)  initiated  by  Iq 

iii)  second  harmonic  ECRM  interaction 

iv)  beam  further  modulated  by 

3rd.  stage:  i)  =  Iq+I^-j  +1^2  + ....  at  the  entrance  and  f^  =  2  f2  =  4  fj 

ii)  Vf4  =  J  E4  “  d(re0)  initiated  by 

iii)  fourth  harmonic  ECRM  interaction 

iv)  Vf 4 ;  I  output  power 
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1 .  Ku  band  fundamental  gyroklystron  2.  Gyroklystron  with  fundamental  prebunch  and 
second  harmonmic  output  3.  Gyroklystron  with  second  haimonic  prebunch  and  fourth 
hannonic  output  4.  fourth  harmonic  W  band  gyroklystron  amplifier  5.  Phase- 
controlled  harmonic-multiplying  inverted  gyrotwyctron  ( Phigtron ). 


A 

I 
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1.  Argon  welding  edge  2.  Cut-off  waveguide  3.  Tuning  plug  assembly  4. 
Input  window  support  5.  Input  window  assembly  6. Ceramic  tube  T.Antenna 
coil  8.  Body  sleeve  9.  Tg<OTEoi<->TEo2  Mode  convertor  10.  NoftTEo2 
mode  suppressor  1 1 .  Mode  filter  1 2.  Attenuator  13.  Transition  waveguide 


6^KG 
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Table  1*  Design  Parameters  of  the  Proof-of-Principle  Experiment 


Phigtron 

Gyro-TWT 

Output  center  frequency 

35  GHz 

35  GHz 

Input  center  frequency 

17.5  GHz 

17.5  GHz 

Bandwidth 

>50  MHz 

>5% 

(phase-locking) 

(instantaneous) 

Output  power 

>  200KW(peak) 

>  100  KW  (peak) 

Efficiency 

>30% 

>15% 

Gain 

40  dB 

30  dB(linear) 

(phase-locking) 

25  dB  (saturated) 

Harmonic  number  s 

2 

2 

Length  of  interaction  region 

40  cm 

50  cm 

Output  mode 

Magnetic  field 

6.5  kG(max.)' 

6.5  kG(max.)‘ 

Gun  type 

MIG' 

MIG' 

Gun  voltage 

60  kV 

60  kV 

Gun  current 

12  A 

12  A 

'  Compatible  with  modem  permanent  magnets  ( Nd/Fe/B). 
*  Magnetron  Injection  Gun. 


■  Based  on  the  cold-test  and  calculation  results  (K.R.Chu  etal.  to  be  published)  we  have  designed  a  two- 
mge  phigtron  oscillator  and  a  two-stage  gyro-TWT  amplifier  for  proof-of-principle  demonstration. 
These  devices  will  soon  be  in  hot-test  experiment. 
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4  Efficiency  Enhancement  of  Gyrotrons  by  Using  De¬ 
pressed  Collectors 

High  power  millimeter  waves  sources,  such  as  Gyrotrons  and  Free  Electron  Lasers  do  not 
presently  have  the  efficiency  of  conventional  tubes  operating  at  longer  wave  lengths.  Over¬ 
all  efficiency  is  an  important  consideration  in  various  applications  such  as  plasma  heating, 
materials  processing,  etc.  The  Institute  of  Plasma  Research  has  been  actively  engaged  in 
studies  on  enhancement  of  overall  efficiency  by  using  depressed  collectors. 

Techniques  have  been  developed  for  tailoring  the  electric  and  magnetic  field  configurations 
in  the  collector  region,  such  that  electrons  in  the  spent  beam  would  be  sorted  according  to 
energy  and  would  not  turn  back  before  getting  collected.  These  techniques  have  been  applied 
to  the  cases  of  large  and  small  orbit  gyrotrons.  Designs  were  developed  for  tubes  with  axially 
and  radially  extracted  beams,  and  for  cavity  type  as  well  as  quasi-optical  type  resonators 
for  rf  interaction. 

A  library  of  computer  codes  has  been  developed  to  assist  in  various  stages  of  depressed 
collector  design.  A  flowchart  is  given  in  Annexure  A  and  a  description  of  the  codes  is  given 
in  Annexure  B.  They  help  to  define  the  regions  of  accessibility  for  electrons,  and  provide  pre 
and  post  processors  for  the  code  which  simulates  electron  trajectories.  They  help  to  reduce 
the  design  time  and  evaluate  new  designs  in  terms  of  collector  efficiency  and  heat  dissipation 
density  profile.  Among  other  things  they  provide  insights  into  the  behavior  of  the  electron 
trajectories  through  three  dimensional  representation  of  an  effective  potential  (Figure  1). 
Animated  contours  can  be  generated  to  get  an  overview  of  the  effect  of  variations  due  to 
different  parameters.  A  summary  of  the  work  was  presented  at  the  International  Conference 
on  Infra  Red  and  Millimeter  Waves  in  Sendai,  Japan,  in  December  1994.^ 

A  feasibility  study  was  done  for  a  single  stage  depressed  collector  for  a  gyrotron  for 
the  International  Thermonuclear  Experimental  Reactor  (ITER).  Currently  a  design  is  being 
developed  for  a  two  stage  depressed  collector  for  it  and  for  a  110  GHz  gyrotron,  in  close 
collaboration  with  industry.  The  library  of  codes  has  been  provided  to  Calabazas  Creek 
Research  in  the  form  of  diskettes  and  user  manuals  for  the  software.  Consultancy  is  being 
provided  to  them  for  the  design  of  two  depressed  collectors  for  a  millimeter  wave  gyrotron 
with  power  output  of  1  Megawatt.  There  is  also  close  collaboration  with  the  tube  division 
of  Communications  and  Power,  Inc.  (Ex  Varians),  for  this  project. 


^Amarjit  Singh,  Daniel  S.  Weile  and  Victor  L.  Granatstein.  “Simulation  and  Comparison  of  Depressed 
Collector  Configurations  for  a  Gyrotron”.  1994  Conference  Record  IR  &  MM  Waves,  pp  73-74,  Sendai, 
Japan. 
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Psi**' 


3-D  Plot  of  Effective  Potential  for 
Two  Depressed  Collectors 


Figure  1 
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Flow  Diagram  for  Design  of  Depressed  Collectors 


No 


Rerun  'POISSON' 
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Yes 

currents 

Annexure  A 
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Computer  Codes  for  Design  of  Depressed  Collectors 

At  the  Institute  for  Plasma  Research,  University  of  Maryland,  we 
have  developed  a  library  of  computer  codes  as  aids  to  design  of 
depressed  collectors,  especially  for  gyrotrons.'  The  tasks  for 
which  we  have  codes  now  are: 

*  !Raypick3*  for  preparation  of  data  arrays  regarding  the 
spent  beam  for  input  to  the  code  for  trajectory  tracing.  The  basic 
data  are  derived  from  the  simulation  of  rf  interaction.  Criteria  can 
be  inserted  for  sampling  and  randomization,  and  calculations 
made  for  the  parameters  needed  by  the  trajectory  code,  'Egun'. 

*  'Spectrum'  for  evaluation  of  the  current  and  power  spectra 
as  a  function  of  energy  for  a)  the  spent  beam  as  characterised  by  rf 
simulation,  and  b)  the  representative  set  of  rays  as  given  by 
'Raypick'. 

*  'ProfilEM*  for:  a)  precise  control  of  magnetic  flux  lines  in 
depressed  collector  region,  b)  evaluation  of  contours  of  effective 
potential  in  the  collector  region,  and  making  2-D  and  3-D  plots  of 
the  same,  so  as  to  ensure  collection  of  all  rays  and  to  give  energy 
sorting  and  c)  pre-processors  for  linking  this  code  to  the 
output/input  of  the  codes,  'Poisson',  and  'Egun' 

*  ITesteff  for  calculation  of  collector  efficiency  from  an 
output  file  in  'Egun'. 

*  TidissS'  for  evaluation  of  the  profile  of  heat  dissipation 
density  on  the  surface  of  the  collectors.  The  strip  width  for  which 
the  density  is  calculated  can  be  varied.  The  heat  dissipation 
density  profiles  are  of  particular  interest  for  elimination  of 
potentially  harmful  hot  spots. 


Annexure  B 
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Features  of  ProfilEM 

ProfllEM  was  programmed  with  a  visual  interface  in 

Windows™.  You  enter  data  and  make  choices  with  user 

friendly  devices  such  as  dialog  boxes,  option  buttons  and  action 

buttons.  Specifically,  ProfllEM  can: 

•  interactively  adjust  contour  lines  and  magnetic  field 
profiles  by  allowing  you  to  change  currents,  potentials  and 
canonical  angular  momentum  and  show  the  results  of  your 
changes  in  a  few  seconds, 

•  display  contours  and  profiles  against  the  backdrop  of  the 
device  geometry  so  that  the  configuration  of  the  fields  vis- 
a-viz  the  geometry  may  be  optimized, 

•  zoom  into  a  narrower  viewing  area  to  emphasize  a  certain 
region  of  interest, 

•  create  hard  copies  of  designs  with  data  on  the  electrode 
potentials,  coil  currents,  and  canonical  angular  momenta 
incorporated  in  the  design  so  that  you  can  keep  a  record  of 
the  design  as  it  evolves, 

•  create  text  files  which  contain  arrays  that  can  be  imported 
by  other  programs  (such  as  Mathematica’’’’^)  to  draw  three 
dimensional  plots  of  the  potential  surface, 

•  make  script  files  so  that  universal  problem  parameters  only 
need  to  be  defined  by  the  user  once  and  subsequently  can 
be  loaded  from  disk, 

•  make  and  show  film  loops  that  show  how  contours  and 
profiles  evolve  with  changing  currents,  potentials  or  the 
canonical  angular  momentum. 
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TASK  B:  Prebunched  Gyrotron  Studies 


1  Overview 

In  recent  years,  numerous  theoretical  and  experimental  investigations  of  gyrotrons  operating 
near  the  fundamental  cyclotron  frequency  have  proven  them  to  be  reliable,  efficient,  high 
power  sources  of  microwave  and  millimeter  wave  radiation  [1,  2].  Unfortunately,  a  serious 
drawback  of  the  classic  gyrotron  configuration  at  high  frequencies  is  the  high  magnetic  field 
level  required  for  first  harmonic  operation  and  the  associated  power  consumption.  Thus, 
they  are  not  viable  candidates  for  some  high  frequency  applications  that  require  compact 
and  lightweight  tubes.  One  way  to  decrease  the  requisite  magnetic  field  is  to  operate  at 
an  harmonic  of  the  cyclotron  frequency.  Unfortunately,  the  efficiencies  of  second  harmonic 
devices  usually  lag  behind  their  first  harmonic  counterparts  by  several  percent  and  efficiencies 
tend  to  decrease  fairly  rapidly  with  increasing  harmonic  number  after  that  point. 

We  have  recently  proposed  a  method  of  pre-bunching  an  electron  beam  to  enhance  the 
efficiency  in  high  harmonic,  large-orbit  devices  which  is  based  on  the  cusp-injection  scheme 
[3,  4].  This  device  (shown  in  Fig.  1),  which  we  call  the  axially-modulated,  cusp-injected, 
large-orbit  gyrotron,  uses  the  axial  bunching  mechanism  of  the  klystron  amplifier  on  an 
annular,  linearly-streaming  beam  to  generate  the  pre-bunched  beam  [5].  The  axial  bunches 
are  created  ballistically  after  the  beam  passes  through  an  input  cavity  which  is  driven  with  a 
TM„i„o  circularly-polarized  mode.  The  axial  bunches  are  converted  into  azimuthal  bunches 
when  the  beam  encounters  a  non-adiabatic  magnetic  transition  (cusp).  If  the  drive  frequency 
is  at  the  mth  harmonic  of  the  cyclotron  frequency,  the  net  effect  is  that  the  beam  will  form 
m  spokes  that  will  rotate  at  the  cyclotron  frequency.  If  an  output  cavity  that  operates  in  a 
TEtoii  mode  is  placed  immediately  after  the  cusp,  efficient  energy  extraction  can  occur.  For 
example,  an  earlier  design  study  [5]  demonstrated  that  an  efficiency  of  40%  was  achievable 
in  X-Band  at  the  fourth  harmonic  via  an  interaction  between  a  100  kV,  25  A  electron  beam 
and  a  simple  TE411  cavity. 

While  this  initial  result  is  quite  promising,  there  are  some  practical  limits  to  the  appli¬ 
cable  range  of  parameters  for  this  device.  For  example,  the  beam  radius  of  the  large-orbit 
configuration  is  equal  to  the  Larmor  radius,  which  is  proportional  to  the  perpendicular  veloc¬ 
ity  for  a  given  cyclotron  frequency.  Thus,  the  beam  radius  in  low  voltage  systems  decreases 
fairly  substantially  with  decreasing  beam  energy.  However,  the  relative  radius  of  the  maxi¬ 
mum  electric  field  in  a  cavity  increases  with  increasing  harmonic  number.  In  fact,  the  electric 
field  at  radii  much  less  than  the  maximum  field  radius  varies  by  a  factor  that  is  proportional 
to  the  power  of  the  harmonic  number.  In  such  systems,  the  achievement  of  the  necessary 
field  strength  for  efficient  interaction  often  leads  to  unrealistic  requirements  on  output  cavity 
quality  factors. 

There  are  at  least  two  possible  ways  to  overcome  these  difficulties  and  design  efficient,  low 
voltage,  harmonic  devices.  One  way  is  to  use  a  vane-resonator  to  couple  various  azimuthal 
modes  together.  The  resultant  cavity  modes  can  support  the  necessary  azimuthal  field  vari¬ 
ation  with  a  much  lower  cutoff  frequency.  Thus,  the  wall  radius  can  be  reduced  considerably 
as  compared  to  the  smooth  wall  case  and  the  interaction  impedance  can  be  significantly 
increased.  The  second  approach  is  to  pre-bunch  a  small-orbit  beam.  Axial  modulation  is 
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still  a  viable  concept  because  a  small  orbit  beam  can  be  generated  by  a  non-adiabatic  mag¬ 
netic  transition  that  does  not  reverse  the  direction  of  the  field.  The  guiding  center  radius 
can  be  adjusted  to  place  the  beam  at  a  location  where  the  electric  field  is  relatively  strong. 
This  non-zero  guiding  center  radius  expands  the  number  of  potential  modes  that  the  beam 
can  interact  with.  While  this  may  enable  us  to  select  a  desirable  output  mode,  it  also  can 
potentially  lead  to  additional  mode  competition  problems.  A  properly  pre-bunched  beam, 
however,  should  help  to  reduce  mode  competition  by  preferentially  selecting  the  desired 
operating  mode.  Efficiencies  for  both  types  of  devices  can  significantly  exceed  those  of  non- 
prebunched  beams.  Yet,  many  of  the  stability  problems  usually  associated  with  gyrotron 
amplifiers  are  not  an  issue  with  this  novel  approach. 

The  goal  of  this  program  is  to  design,  build,  and  test  a  proof-of-principle  experiment 
based  on  this  prebunched  gyrotron  concept.  To  be  of  interest  to  applications  that  require 
compact  devices,  we  hmit  our  designs  to  voltages  below  50  kV.  To  demonstrate  moderate 
power  levels,  we  have  focused  on  designs  that  produce  peak  powers  on  the  order  of  100 
kW.  We  have  kept  our  design  frequencies  in  the  5-10  GHz  range  in  order  to  perform  a  cost 
effective  experiment.  The  designs  are  fairly  conservative  and  can  no  doubt  be  pushed  to 
higher  frequencies  with  a  power  that  scales  better  than  the  normal  expectation. 

The  first  year’s  effort  has  been  spent  predominantly  working  on  three  areas.  The  first 
-  effort  was  to  develop  the  necessary  modeling  tools  to  analyze  the  microwave  circuits  which 
can  efficiently  interact  with  low-voltage  rotating  beams  in  a  variety  of  configurations.  The 
five  main  codes  developed  or  modified  for  this  effort  are  described  in  the  following  section. 
The  second  task,  which  represented  the  bulk  of  the  work,  was  to  use  these  new  codes  along 
with  some  existing  codes  to  design  a  proof-of-principle  experiment  that  can  be  carried  out 
in  our  laboratory.  Included  in  the  appendix  are  copies  of  two  papers  [6,  7]  that  have  been 
submitted  to  journals  and  which  describe  the  results  of  our  microwave  circuit  designs.  These 
papers  are  very  briefly  summarized  in  the  next  section,  along  with  descriptions  of  the  design 
work  not  included  in  the  journal  articles.  The  final  task  was  to  begin  the  construction  of 
the  experiment.  The  initial  phase  involved  bringing  additional  personnel  onto  the  project 
(engineers,  students,  etc.),  identifying  existing  hardware  that  was  available  and  compatible 
with  the  designs,  and  procuring  additional  hardware.  These  efforts  are  described  in  the  third 
section. 

In  the  final  section  we  describe  our  proposed  effort  for  the  second  year. 


2  Theoretical  Effort 

2.1  Code  development 

Five  main  codes  were  written  or  enhanced  in  the  past  year  to  assist  in  the  design  of  low 
voltage,  prebunched  gyrotrons  systems.  Two  “driver”  codes  focused  on  the  microwave  input 
system.  The  first  computes  the  field  profiles,  quality  factors,  and  required  drive  power  for 
TMmnp  dielectric-loaded,  right-circular  cavities.  Up  to  three  ideal  dielectric  regions  can  be 
introduced  between  arbitrary  radial  locations.  The  second  code  computes  the  field  profiles  for 
a  sheath-helix  structure.  The  input  system,  which  will  be  used  in  the  design  of  a  broadband 
version  of  this  device,  can  have  a  conducting  wall  inside  or  outside  the  helix  (or  both). 
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The  “buncher”  code  uses  the  output  from  the  “driver”  codes  to  simulate  the  velocity 
modulation  and  ballistic  bunching  of  the  beam.  It  can  also  be  used  without  the  driver 
codes  if  the  input  system  is  a  simple  right-circular  (unloaded)  cavity.  Initial  uniform  particle 
distributions  in  time,  azimuthal  angle,  and  axial  location  can  be  specified.  Finite  beam 
thicknesses  are  modeled  with  three  layers  of  electrons.  A  predictor-corrector  scheme  is  used 
to  numerically  integrate  the  particles  through  an  acceleration  region,  the  input  cavity,  the 
drift  region,  and  the  nonadiabatic  magnetic  transition.  The  transition  uses  a  piecewise- 
linear  model  and  unbalanced  transitions  can  be  specified.  Different  integration  step  sizes 
are  possible  in  the  cavity,  drift,  and  transition  regions.  A  Runga-Kutta  algorithm  is  used 
to  start  the  integration  whenever  necessary.  The  code  outputs  the  phase  space  distribution 
of  the  particles  and  various  statistics.  In  addition  to  calculating  the  average  guiding  center 
position,  average  velocity  ratio,  and  the  perpendicular  and  parallel  velocity  spreads,  the  code 
also  computes  the  RMS  spread  in  beam  energy  and  a  measure  of  the  bunching  effectiveness. 
This  measure  first  computes  the  RMS  spread  in  angular  location,  but  then  normalizes  it  so 
that  a  completely  unbunched  beam  would  yield  a  value  of  100%.  We  normally  report  100% 
minus  this  value,  so  that  an  unbunched  beam  has  a  bunching  effectiveness  of  0%  and  a  fully 
bunched  beam  has  an  effectiveness  of  100%.  Because  particles  can  be  launched  at  various 
times,  the  angular  location  is  given  relative  to  the  phase  of  an  EM  wave  at  the  frequency  of 
-  the  buncher  cavity. 

The  “VR  cavity”  code  computes  the  resonant  frequency,  quality  factor,  and  field  profiles 
of  a  closed  cavity  with  vane  resonators.  The  number  of  terms  kept  in  the  infinite  series 
can  be  specified.  The  usual  approximation  of  a  constant  field  in  the  slot  is  used  to  find  the 
dispersion  relation.  An  arbitrary  phase  shift  between  vanes  can  be  introduced  provided  that 
it  is  consistent  with  the  boundary  conditions.  The  quality  factor  is  computed  by  assuming 
a  good  copper  cavity  and  numerically  integrating  the  surface  currents  on  all  walls. 

The  “efficiency”  code  uses  the  phase-space  results  from  the  buncher  code  to  compute 
the  net  interaction  in  the  output  cavity.  There  are  now  three  ways  to  specify  an  EM  field 
configuration.  Simple  right-circular  field  profiles  can  be  generated  internally  and  the  fields 
of  a  VR  structure  can  be  read  in  from  the  output  of  the  previously  described  code.  Finally, 
we  have  modified  the  code  to  accept  input  from  a  scattering  mati'ix  code  [8].  In  this  way, 
we  can  examine  the  effects  of  non-sinusoidal  field  profiles  and  post-cavity  beam-microwave 
interactions  on  tube  performance.  The  efficiency  code  uses  the  same  numerical  methods 
as  the  buncher  code.  A  linearly-tapered  magnetic  field  can  be  introduced  into  the  region. 
Efficiency  is  computed  as  a  function  of  the  EM  wave’s  amplitude  and  phase.  These  two 
quantities  can  be  swept  automatically  to  aid  in  the  search  for  optimal  efficiency.  The  code 
calculates  the  output  cavity  quality  factor  required  to  achieve  the  specified  fields. 


2.2  Small-orbit  design 

The  results  of  our  small-orbit  design  are  discussed  in  great  detail  in  the  first  paper  contained 
in  the  Appendix.  In  this  section  we  highlight  only  the  key  points  of  the  design. 

The  parameters  of  the  second  harmonic,  small-orbit  design  are  given  in  center  column 
of  Table  1.  A  picture  of  the  system  with  a  typical  beam  profile  is  given  in  Fig.  2.  The 
smooth- walled  output  cavity  has  a  main  section  length  of  10  cm  and  a  coupling  lip  with 
a  radius  of  1.7  cm.  Simulations  are  carried  out  with  the  codes  described  in  the  previous 
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Table  1;  Amplifier  system  parameters. 


Beam  Parameters 

Small-orbit 

Large  orbit 

Voltage  (kV) 

35 

45 

Current  (A) 

4 

4 

Average  radius  (cm) 

1.18 

1.13 

Beam  thickness  (mm) 

7.33 

0.8 

Velocity  ratio  (ux/^^z) 

2.28 

2.20 

Input  Cavity 

Drive  frequency  (GHz) 

9.9 

5.5 

Operating  mode 

TMo20 

TM310 

Resistive  Quality  factor  (Q) 

2800 

2170 

Radius  (cm) 

2.66 

1.701 

Length  (cm) 

0.20 

0.23 

Magnetic  Field  Parameters 

Buncher  guide  field  (kG) 

0.985 

-0.57 

Output  magnetic  field  (kG) 

1.83 

0.61 

Cyclotron  frequency  (GHz) 

4.794 

1.577 

Cusp  width  (cm) 

0.50 

0.38 

Cavity  -  cusp  spacing  (cm) 

13 

30 

Output  field  taper 

-3% 

13% 

Output  cavity 

Output  frequency  (GHz) 

9.9 

5.5 

Operating  mode 

TEoii 

TT-mode 

Diffractive  Q 

1825 

1240 

Resistive  Q 

28,580 

13,430 

Radius  (cm) 

1  1.864 

1.763 

Length  (cm) 

11.09 

15.00 

Harmonic  number 

2 

3 

section  and  are  checked  with  a  modified  gyroklystron  code  [9].  With  a  beam  thickness  of 
6.7%  (prior  to  the  magnetic  transition),  the  net  axial  velocity  spread  at  the  output  cavity 
entrance  is  12%  and  the  net  energy  spread  is  3.2%.  A  slight  decrease  in  the  magnetic  field 
strength  throughout  the  output  region  is  required  to  maximize  efficiency.  The  required  drive 
power  is  840  W  and  an  output  power  of  over  74  kW  is  achieved  with  an  efficiency  of  53%. 
Hence,  the  large-signal  gain  is  approximately  19.5  dB.  The  system  drive  curve  is  shown  in 
Fig.  3.  Powers  as  high  as  78  kW  can  be  achieved  with  drive  powers  of  about  1.2  kW.  The 
zero-drive  efficiency  is  less  than  2%.  An  examination  of  the  start-oscillation  currents  for  all 
the  potential  spurious  modes  indicates  that  only  the  desired  operating  mode  is  unstable  at 
the  nominal  operating  point.  The  only  other  mode  to  have  a  start  current  below  10  A  is  the 
TEin  mode,  but  this  is  zero-drive  stable  for  all  magnetic  fields  above  1.79  kG.  Simulations 
indicate  that  efficiency  is  fairly  insensitive  to  beam  thickness.  They  also  predict  efficiency 
enhancement  of  a  few  percent  for  larger  average  velocity  ratios  (up  to  a  =  3). 
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2.3  Large-Orbit  Design 

Again,  the  results  of  our  large-orbit  design  aie  discussed  in  considerable  detail  in  the  second 
paper  supplied  in  the  Appendix.  In  this  section  we  summarize  the  design  parameters  at  the 
nominal  operating  point  and  discuss  operation  at  other  harmonics. 

The  design  parameters  of  the  third  harmonic,  large-orbit  design  are  given  in  the  final 
column  of  Table  1.  The  output  cavity  is  of  the  vane-resonator  type  and  has  six  slots,  each 
with  a  depth  of  0.75  cm  and  a  relative  opening  of  50%.  The  cavity  cross-section  is  plotted 
in  Fig.  4  along  with  beam  distributions  at  the  entrance  (Fig.  4a)  and  exit  (Fig.  4b)  of  the 
output  cavity.  The  input  cavity  is  assumed  to  contain  lossless  dielectric  material  to  enhance 
the  axial  electric  field  at  the  beam.  With  a  beam  thickness  of  6.7%,  the  net  axial  velocity 
spread  at  the  output  cavity  entrance  is  9.0%  and  the  net  energy  spread  is  3.8%.  A  fairly 
substantial  increase  in  the  magnetic  field  strength  throughout  the  output  region  is  required 
to  maximize  efficiency.  The  required  drive  power  is  1.1  kW  and  an  output  power  of  nearly 
100  kW  is  achieved  with  an  efficiency  of  55%  (and  a  resultant  large-signal  gain  of  nearly  19.6 
dB).  Simulations  indicate  that  this  device  is  considerably  more  sensitive  to  beam  thickness 
than  the  small-orbit  device.  On  the  other  hand,  the  post-cusp  velocity  spread  induced  by 
the  buncher  cavity  is  considerably  less  for  the  large-orbit  beam  than  for  the  previous  system. 

We  have  also  taken  a  preliminary  look  at  efficiency  enhancement  in  higher  harmonic 
systems  with  comparable  beam  voltages.  In  Fig.  5a,  we  plot  the  simulated  efficiency  of 
a  large-orbit  device  as  a  function  of  harmonic  number.  The  point  on  the  dashed  line  at 
the  third  harmonic  corresponds  to  the  system  described  in  the  previous  section.  For  the 
other  harmonic  results,  the  magnetic  field  and  beam  parameters  are  fixed  at  the  values 
given  in  Table  1.  However,  the  cavity-cusp  spacing,  and  the  cavities’  harmonic  number, 
radial  dimensions,  and  number  of  (output)  vanes  are  adjusted  to  the  appropriate  values. 
Furthermore,  the  output  cavity  quality  factor  is  adjusted  to  optimize  efficiency  (Fig.  5b). 
The  figure  also  indicates  the  simulated  efficiency  achievable  for  a  zero  thickness  beam.  With 
the  exception  of  the  second  harmonic  result,  the  finite  beam  thickness  results  in  an  efficiency 
decrease  of  about  10  -  12%. 

Figure  5  clearly  demonstrates  the  rapid  decline  in  tube  efficiency  with  harmonic  number 
which  is  inherent  in  these  systems.  This  problem  is  exacerbated  by  the  dramatic  rise  in  the 
required  quality  factor.  For  example,  the  sixth  harmonic  design  presented  here  requires  a 
quality  factor  over  three  times  the  resistive  Q  of  a  copper  cavity.  While  a  full  optimization  of 
system  design  at  each  harmonic  (e.g.  adjustment  of  the  output  cavity  length,  magnetic  field 
tapering,  etc.)  would  undoubtedly  improve  the  indicated  performance  somewhat,  the  fourth 
harmonic  is  probably  the  highest  viable  harmonic  number  for  the  45  kV  beam  considered 
here. 


2.4  Broadband  bunching  investigation 

An  extension  of  this  device  that  we  have  just  started  to  examine  is  designed  to  produce  similar 
efficiency  enhancement  in  a  broadband  configuration.  The  basic  idea  is  to  replace  the  input 
cavity  with  a  slow  wave  structure  and  the  output  cavity  with  a  gyrotwystron-like  output 
waveguide  [10].  Initial  results  with  the  beam  parameters  of  the  small-orbit  configuration 
indicate  that  bunching  over  a  wide  range  of  frequencies  is  possible  with  only  a  slight 
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Table  2:  The  electron  gun  parameters. 


Parameter 

Value 

Average  cathode  radius  (cm) 

3.0 

Cathode  thickness  (mm) 

2.0 

Cathode  loading  (A/cm^) 

1.0 

Cathode  magnetic  field  (G) 

260 

Magnetic  compression 

4.0 

Peak  electric  field  (kV/cm) 

24 

increase  in  the  resultant  velocity  spread  as  compared  to  the  input  cavity  approach.  Sheath 
helices  can  also  be  used  in  principle  to  excite  the  fields  required  for  the  large-orbit  configura¬ 
tion.  We  are  now  in  the  process  of  modifying  our  codes  to  model  the  traveling-wave  output 
structure. 

2.5  Electron  gun  design 

We  have  been  considering  three  slightly  different  electron  gun  configurations  in  the  course 
of  this  investigation.  EG  UN  [11]  was  used  for  all  simulations.  The  first,  pictured  in  Fig.  6, 
is  a  gun  designed  to  produced  2  A  at  35  kV.  The  second  configuration  is  similar  in  shape 
but  has  a  decreased  anode-cathode  gap  and  is  designed  to  operate  at  4  A  and  45  kV.  The 
space-charge  limiting  current  is  simulated  to  be  6  A.  A  summary  of  the  gun  design  values 
for  the  first  configuration  is  given  in  Table  2.  The  peak  electric  field  on  the  cathode  is  24 
kV /cm. 

We  are  currently  working  on  a  higher  compression  design  that  can  be  used  to  power  the 
third  harmonic  large-orbit  design.  It  is  believed  that  the  same  cathode  can  be  used  for  both 
experiments. 

3  Experimental  Effort 

With  the  completion  of  the  two  design  studies  presented  in  the  previous  section,  we  have 
just  recently  turned  our  attention  to  the  proof-of- principle  experiment.  Our  initial  effort 
to  assemble  a  test  bed  has  included  work  on  a  modulator,  the  magnetic  field  system,  the 
vacuum  system,  microwave  components,  and  test  equipment.  Progress  in  each  of  these  areas 
is  described  in  the  following  paragraphs. 

We  have  obtained  a  PFN-type  modulator  that  has  a  maximum  peak  power  of  about  950 
kW  and  a  pulse  duration  of  either  0.5  fis  or  5  fis.  We  will  use  the  modulator  in  the  long 
pulse-length  mode.  We  have  repaired  some  minor  problems  and  characterized  the  modulator 
performance  as  a  function  of  load  impedance.  Optimal  performance  is  achieved  with  a  pulse 
voltage  of  49  kV  and  a  current  of  about  19  A  (at  90%  of  the  supply’s  maximum  output). 
Since  the  beam  current  is  in  the  2-6  A  range,  a  shunt  resistor  will  be  used  to  match  the  PFN 
impedance. 

We  have  identified  over  a  half-dozen  magnetic  field  coils  that  are  available  to  us  and  that 
will  be  easily  able  to  handle  the  field  requirements  of  the  large-orbit  design.  We  are  just 
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starting  to  evaluate  their  ability  to  work  at  the  higher  magnetic  field  required  for  the  small- 
orbit  design.  If  they  prove  to  be  inadequate,  we  will  attempt  to  procure  several  water-cooled 
magnets.  We  have  also  identified  a  few  power  supplies  that  can  be  used  to  energize  some  of 
the  coils  but  have  ordered  three  more  supplies.  We  have  obtained  an  iron  plate  that,  with  a 
few  modifications,  should  prove  adequate  for  the  non-adiabatic  transitions  of  both  designs. 

We  have  obtained  one  60  £/s  ion  pump/controller  set  and  have  access  to  roughing  and 
bake-out  pumps  in  addition  to  leak-checking  hardware.  We  are  in  the  process  of  ordering 
additional  vacuum  hardware. 

We  have  identified  sources  that  are  available  to  us  and  cover  the  range  from  8.5-10  GHz. 
Two  of  the  sources  are  pulsed  magnetrons  with  considerably  more  power  than  required  but 
pulse  lengths  of  only  about  2  /rs.  Another  source  (near  9.9  GHz)  produces  a  few  kilowatts 
and  can  be  run  CW.  We  intend  to  order  a  5.4  -  5.8  GHz  magnetron  that  can  be  integrated 
into  an  existing  power  supply  and  used  for  the  large-orbit  design.  We  are  also  in  the  process 
of  ordering  a  selection  of  microwave  components  that  will  be  required  to  perform  the  proof- 
of- principle  experiment. 

Finally,  we  have  begun  to  accumulate  the  necessary  test  equipment  that  will  be  required 
to  successfully  complete  this  study  (e.g.  a  computer,  interface  hardware,  oscilloscopes,  etc.) 

4  Proposed  Work  for  the  Next  Year 

We  have  developed  an  ambitious  schedule  for  next  year  which  should  allow  us  to  begin 
microwave  amplifier  experiments  by  the  beginning  of  the  third  year.  We  intend  to  first  try 
to  build  and  test  the  second  harmonic  small-orbit  design,  primarily  because  more  of  the 
hardware  for  this  design  is  already  on  hand  in  our  Institute.  With  the  exception  of  the 
required  circuit  magnetic  field,  the  small-orbit  components  are  somewhat  easier  to  construct 
than  the  corresponding  large-orbit  components.  We  will  continue  to  improve  our  modeling 
capabilities  and  explore  other  pre-bunched  gyrotron  designs.  We  have  briefly  summarized 
the  major  tasks  proposed  for  next  year  below. 

1.  Modify  the  “efficiency”  code  to  automatically  search  for  the  maximum  efficiency  as  a 
function  of  field  amplitude  and  phase. 

2.  Modify  the  “buncher”  code  to  accept  input  from  the  scattering  matrix  code  to  allow 
coaxial  input  cavities  and  to  investigate  the  effects  of  fringing  fields. 

3.  Modify  the  “output”  code  to  model  traveling  wave  interactions. 

4.  Modify  the  “helix”  and  “buncher”  codes  for  non-symmetric  modes. 

5.  Modify  the  “buncher”  code  to  accept  input  from  EGUN  simulations. 

6.  Use  HFSS  to  model  the  input  cavities. 

7.  Design  and  evaluate  broadband  amplifier  configurations. 

8.  Re-evaluate  previous  designs  with  the  improved  codes. 
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9. 

10. 

11. 

12. 

13. 

14. 

15. 

16. 

17. 

18. 


Finalize  electron  gun  design,  order  emitter  strip,  and  construct  remaining  components. 
Complete  procurement  of  experimental  hardware. 

Construct  and  cold-test  input  and  output  cavities  for  small-orbit  design. 

Design  and  construct  vacuum  jacket.  Utilize  a  gate-valve  if  at  all  possible  to  keep  gun 
under  vacuum  during  tube  rebuilds. 

Construct  complete  magnetic  field  system  and  characterize  field  profile  in  all  regions, 
particularly  near  the  iron  plate. 

Construct  the  vacuum-compatible  microwave  circuit. 

Construct  beam  diagnostics. 

Characterize  beam  performance. 

Construct  microwave  diagnostics. 

Initiate  microwave  amplifier  experiment  (most  likely  not  to  begin  until  the  third  year). 
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TASK  C: 

Studies  of  Cyclotron  Maser  Amplifiers  with 
Trochoidal  Electron  Beams  (Trochotrons) 


After  developing  the  linear  theory  of  trochotrons  the  quality  of  trochoidal  electron  beams 
has  been  studied.  Two  concrete  designs  of  electron  guns  have  been  analyzed  intended  for 
operation  at  35  GHz  in  second  harmonic  trochotrons.  It  was  shown  that  trochoidal  electron 
beams  can  be  formed  with  large  ratios  of  perpendicular  to  parallel  (drift)  velocity  and  very 
small  spread  in  drift  velocities.  This  makes  such  beams  suitable  for  utilization  in  trochotron- 
traveling-wave  amplifiers. 

A  paper  submitted  for  publication  in  the  6th  Special  Issue  of  the  IEEE  Transactions  on 
Plasma  Science  on  High-Power  Microwaves  is  included  as  an  Appendix. 
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APPENDIX 


TASK  A: 


TASK  B: 


TASK  C: 


Papers  Pertaining  to  this  Research 


“Phase-locking  of  a  second-harmonic  gyrotron  oscillator 
using  a  quasi-optical  circulator  to  separate  injection  and 
output  signals” 


“Design  of  a  high  efficiency,  low  voltage,  axially  modulated, 
cusp-injected,  second  harmonic,  X-band  gyrotron  amplifier” 

“Design  of  an  efficient,  low  voltage,  third  harmonic,  large-orbit 
gyrotron  amplifier  with  a  vane-resonator  output  cavity” 


“On  the  design  of  electron  guns  for  cyclotron  resonance  masers 
utilizing  trochoidal  electron  beams  (trochotrons)” 
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Phase-Locking  of  a  Second-Harmonic  Gyrotron  Oscillator 
Using  a  Quasi- Optical  Circulator  to  Separate  Injection  and 

Output  Signals* 

H.Z.  Guo,  D.J.  Hoppe,  Member,  IEEE,  J.  Rodgers,  R.M.  Perez,  J.P.  Tate,  Member,  IEEE, 
B.L.  Conroy,  V.L.  Granatstein,  FELLOW,  IEEE,  A.  Bhanji,  P.E.  Latham,  Member,  IEEE,  ’ 
G.S.  Nusinovich,  Senior  Member,  IEEE,  M.  Naiman,  and  S.H.  Chen 


Abstract — Phase-locking  in  a  34.5  GHz  special  complex  cav¬ 
ity  gyrotron  oscillator  operating  at  the  second  harmonic  of 
the  electron  cyclotron  frequency  was  studied.  Injection  of 
the  locking  power  was  made  via  a  quasi-optical  circulator 
connected  to  the  gyrotron  output.  locking  bandwidth  was 
measured  by  comparing  the  phase  of  the  injection  signal 
and  output  signal  using  a  balanced  mixer.  Locking  was  ob¬ 
served  with  input  power  level  as  low  as  40  dB  below  the  gy¬ 
rotron  output  power.  The  locking  bandwidth  is,  however, 
narrower  than  in  gyrotrons  operating  at  the  fundamental 
cyclotron  frequency  which  may  be  attributed  to  the  longer 
resonant  cavity  in  the  second  harmonic  gyrotron  and  the 
corresponding  larger  value  of  external  quality  factor.  The 
measurements  are  roughly  in  agreement  with  predictions  of 
Adler’s  phase-locking  equation  which  is  given  for  our  system 
in  terms  of  powers  propagating  in  the  output  waveguide  to¬ 
ward  and  away  from  the  gyrotron  cavity. 


I.  Introduction 


Efficient,  high  power,  coherent  millimeter  wave  sources 
are  currently  needed  for  satellite  and  space  communications 
and  for  imaging  radar  to  be  used  in  space  exploration  [1]. 
A  weakly  relativistic  (electron  energy  <  100  keV)  gyrotron 
phase-locked  oscillator  operating  at  high  harmonic  of  the 
electron  cyclotron  frequency  (harmonic  number  5  >  2)  has 
been  considered  as  one  of  the  appropriate  source  choices 
for  these  applications  for  the  following  reasons: 

(1)  the  transverse  dimensions  of  both  the  electron  beam 
and  the  circuit  can  be  much  larger  than  in  klystrons 
or  slow-wave  microwave  tubes  because  of  the  ability  of 
gyrotrons  to  operate  in  fast  wave,  higher-order,  trans¬ 
verse  electric  modes,  resulting  in  higher  average  and 
peak  power  rating  at  a  given  operating  frequency; 

(2)  the  moderate  voltage  requirement  avoids  the  need  for 
large  insulators  or  electron  accelerator  sections;  and 

(3)  the  required  magnetic  field  can  be  reduced  by  the  har¬ 
monic  factors  compared  with  gyrotrons  operating  at 
fundamental  cyclotron  frequency,  potentially  making 
gyro-devices  compatible  with  modern  permanent  mag¬ 
net  technology  up  to  operating  frequencies  of  approx¬ 
imately  (s  X  20)  GHz. 


The  research  described  in  this  paper  wets  carried  out  in  p2Lrt  by  the 
Jet  Propulsion  Laboratory,  California  Institute  of  Technology,  under  a 
contract  with  the  National  Aeronautics  and  Space  Administration;  in  ad¬ 
dition,  this  work  was  supported  in  part  by  the  DoD  Vacuum  Electronics 
Initiative  and  managed  by  the  Air  Force  Office  of  Scientific  Research  un¬ 
der  Grant  AFOSR-91-0390. 
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Recently,  major  advances  have  been  reported  in  gyrotrons 
operating  at  harmonics  of  the  cyclotron  frequency.  High 
efficiency  and  high  power  have  been  achieved  in  free  run¬ 
ning  second  harmonic  gyrotron  oscillators  operating  in  the 
millimeter-wave  band  [2-4].  Furthermore,  there  have  been 
successful  demonstrations  [5,6]  of  phase-locked  gyrotrons 
operating  at  the  fundamental  cyclotron  frequency.  How¬ 
ever,  no  phase-locked  harmonic  gyrotron  operation  had 
been  reported  prior  to  the  present  paper. 

In  this  paper  we  report  the  operation  and  measurement 
of  a  high  power  (~100  kW)  second  harmonic  gyrotron 
phase-locked  by  signal  injection  through  the  output  win¬ 
dow  using  a  quasi-optical  circulator  system  [7]  to  separate 
the  input  and  output  signals.  Section  II  gives  a  brief  de¬ 
scription  of  the  second  harmonic  gyrotron  and  the  arrange¬ 
ment  of  the  phase-locking  experiment.  Section  III  discusses 
the  pertinent  theory  of  phase-locking  taking  into  account 
reflection  of  the  injected  signal  at  the  entrance  to  the  gy¬ 
rotron  cavity;  requirements  on  the  flatness  of  the  voltage 
pulse  applied  to  the  gyrotron  are  also  discussed.  In  Sec¬ 
tion  IV  the  results  of  the  phase-locking  experiment  are  pre¬ 
sented  and  compared  with  theoretical  predictions.  Finally, 
a  summary  is  presented  in  Section  V. 


II.  Description  of  the  Gyrotron  and 
THE  Experimental  Arrangement 


A.  The  Second-Harmonic  Gyrotron 

A  second-harmonic  gyrotron  which  features  a  novel  com¬ 
plex  cavity  as  its  beam-wave  interaction  structure  [2]  was 
used  for  the  demonstration  of  phase-locking.  The  parame¬ 
ters  of  this  gyrotron,  GY-32,  are  shown  in  Table  1.  Figure 
I  is  a  cross-sectioned  drawing  of  the  gyrotron  tube  along 
with  a  typical  magnetic  field  profile.  The  special  complex 
cavity  consists  of  five  parts  as  depicted  in  Fig.  2.  Part  I  and 
part  2  are  uniform  circular  waveguide  sections  with  their 
inner  radii  denoted  by  a  and  6,  respectively.  These  sections 
are  designed  to  operate,  respectively,  in  the  TEqi  and  TE03 
modes.  In  order  to  couple  these  two  modes  into  one  normal 
mode  of  a  complex  cavity  the  eigenfrequencies  of  these  two 
partial  modes  must  be  the  same.  This  is  accomplished  by 
choosing  radii  a  and  6  to  satisfy  the  condition 


P02  fl03  , 


(1) 


w^  ere  eigennumbers  /io2  J>nd  //03  are  roots  of  the  equations 
(/^02)  =  —  0,  Jq  is  the  zero  order  Bessel  func¬ 

tion  and  kc  is  the  cutoff  wavenumber,  which  is  the  same 


for  both  modes.  Then,  if  the  lengths  of  both  sections  are 
also  equal  that  gives  closely  matched  axial  wavenumbers 
for  both  modes;  thus,  one  can  suppose  that  in  a  complex 
cavity  these  two  partial  modes  will  have  closely  matched 
eigenfrequencies. 


Fig.  1.  Cross-sectional  drawing  of  the  GY-32  gyrotron 
tube  along  with  a  typical  magnetic  field  profile:  1)  cath¬ 
ode,  2)  control  anode,  3)  second  anode,  4)  complex  cavity, 
5)  beam  collector,  6)  water  cooling  jacket,  7)  output  win¬ 
dow,  8)  main  magnetic  field  coils,  9)  gun  coils. 
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Fig.  2.  Special  complex  cavity  with  internal  mode  con¬ 
verter  for  transforming  TE02  to  TE03. 

To  provide  a  coupling  between  these  two  modes  a  slot¬ 
ted  waveguide  with  a  gradually  changing  cross  section  is 
used  (part  3).  This  part  is  composed  of  a  main  waveguide 
of  radius  a  and  N  slotted  branch  waveguides.  The  angle 
subtended  by  each  slot,  20o,  varies  slowly  from  0  to  2t/A^ 
over  the  length  cd]  the  depth  of  each  branch  waveguide  is 
constant  and  equal  to  (6  —  a).  Thus,  the  gradually  chang¬ 
ing  slotted  waveguide  efficiently  converts  a  TE02  mode  to 


a  TEo3  mode.  The  three  parts  described  above  together 
with  part  4  (cutoff  section)  and  part  5  (diffraction  output 
section)  form  a  resonator.  Good  mode  selectivity  results 
from  this  kind  of  complex  cavity  for  the  following  reasons: 

(1)  The  radial  metallic  fins  at  one  end  of  the  mode  con¬ 
verter  function  as  a  mode  filter.  This  is  because  modes 
which  have  a  radial  component  of  rf  electric  field  may 
not  be  excited  in  the  presence  of  the  radial  metallic 
sheets. 

(2)  The  mode  converter  functions  as  mode  isolator.  With 
the  design  criterion  applied  to  this  case,  only  one  pair 
of  modes  (TE02  and  TE03)  will  be  coupled  and  form 
a  normal  mode  for  the  whole  system. 

Table  1.  Performance  parameters  of  the  GY-32  gyrotron. 


Harmonic  Number 

2 

Beam  Voltage 

40-60  kV 

Frequency 

34.5  GHz 

Beam  Current 

8-11  A 

Peak  Output  Power 

80-230  kW 

CC  =  t;_L/v|j 

1.5 

Efficiency 

~35% 

Quality  Factor,  Q 

2000 

Ohmic  Quality  Factor,  Qohm 

10,000 

ModeTurity  (TE03  mode) 

>90% 

Due  to  the  superior  mode  selectivity  of  the  special  com¬ 
plex  cavity,  stable,  efficient,  harmonic  gyrotron  operation 
and  high  output  mode  purity  result.  Both  of  these  are  vital 
for  realization  of  the  proposed  phase-locking.  A  calibrated 
TEon  water  calorimeter  placed  at  the  output  of  the  gy¬ 
rotron  verified  that  greater  than  95%  (limit  of  calorimeter 
accuracy)  of  the  power  was  in  the  symmetric  TE  modes. 
Spectral  analysis  of  the  output  power  also  showed  that 
the  TEo3  signal  was  at  least  25  dB  above  any  other  TE 
modes  resulting  in  an  estimated  TE03  purity  of  greater 
than  90%.  It  would  have  been  more  difficult  to  achieve  the 
required  stability  of  operation  in  a  more  usual  smooth-wall 
cavity  gyrotron  operating  at  a  harmonic  of  the  cyclotron 
frequency  because  of  mode  competition  problems.  Also, 
the  phase-locking  arrangement  in  the  present  study  em¬ 
ployed  a  quasi-optical  circulator  for  separation  of  injection 
signal  from  the  gyrotron  output  signal,  and  this  circulator 
has  critical  demands  on  gyrotron  output  mode  purity  for 
its  effective  operation. 

Evidence  for  the  stability  of  gyrotron  operation  is  shown 
in  Fig.  3.  Figure  3(a)  shows  a  typical  modulator  voltage 
and  current  pulse;  3(b)  is  a  detected  sample  of  the  gyrotron 
rf  output;  3(c)  is  the  gyrotron  output  frequency  spectrum; 
and  3(d)  shows  the  output  mode  pattern  recorded  on  a 
liquid  crystal  display.  All  of  the  above  traces  were  taken 
prior  to  injecting  the  phase-locking  signal  (i.e.,  the  gyrotron 
was  in  a  free-running  condition).  Since  gyrotron  shot-to- 
shot  stability  and  instantaneous  frequency  stability  were 
of  great  importance,  gyrotron  output  power  was  traded  off 
for  increased  stability  by  keeping  the  beam  current  and  the 


beam  voltage  at  the  lower  end  of  the  values  shown  in  Ta¬ 
ble  1.  Lower  current  reduces  beam  space  charge  thereby 
increasing  the  coherence  of  the  electron  cyclotron  phase 
bunching.  Lowering  the  voltage  had  the  effect  of  reducing 
shot-to-shot  pulse  height  fluctuations  in  the  modulator  to 
approximately  one  percent.  The  gyrotron  system  operat¬ 
ing  parameters  for  the  phase-locked  experiment  are  shown 
in  Table  2. 


Fig.  3.  Stable  operation  of  the  GY-32  gyrotron.  (a) 
Top  trace  is  the  modulator  voltage,  10  kV/div;  bottom 
trace  is  the  beam  current,  1  A/div;  horizontal  time  axis, 
1  psec/div.  (b)  Detected  gyrotron  output  pulse,  (c)  Spec¬ 
trum  analyzer  display  of  gyrotron  output  from  IIP11970A 
waveguide  mixer.  The  sidebands  shown  are  the  second  ax¬ 
ial  mode  of  the  complex  cavity  generated  during  the  rise  of 
the  beam  pulse,  (d)  Mode  pattern  produced  by  gyrotron 
output  in  TEo3  mode. 


Table  2.  Operating  parameters  for  the  phase-locked 
gyrotron  experiment. 


Beam  Voltage 

42-45  kV 

Beam  Current 

8  Amps 

Frequency 

34.52  GHz 

Output  Power 

80-100  kW 

Gun  Field 

1.15  kG 

Cavity  Field 

6.5  kG 

Repetition  Rate 

43  Hz 

Pulse  Width 

7  /is 

B.  Experimental  Arrangement 

A  block  diagram  of  the  phcise-locked  system  is  shown  in 
Fig.  4.  A  synthesized  microwave  generator  amplified  by  a 
100-watt  traveling  wave  tube  amplifier  (TWTA)  provided 
the  injection  signal.  The  signal  transmission  path  included 
a  conventional  waveguide  isolator  providing  more  than  35 
dB  of  isolation,  a  directional  coupler  for  sampling  the  injec¬ 
tion  signal,  and  a  smooth-wall  dual-mode  horn  providing  a 
gaussian  beam  for  launching  into  the  cjuasi-optic  system. 

The  quasi-optic  circulator  is  depicted  in  detail  in  Fig.  5. 
The  microwave  beam  was  reflected  off  a  flat  plate,  fo¬ 
cused  by  an  ellipsoidal  mirror  with  focal  lengths  of  24  and 
48  inches,  and  then  redirected  towards  a  polarized  grid. 
The  grid  was.  fabricated  by  printing  0.00615  inch  copper 
strips  on  0.0123  inch  centers  on  a  0.189  inch  thick  disk  of 
fused  quartz.  The  beam  then  passed  through  the  circula¬ 
tor  which  consists  of  a  5.75  inch  diameter  disk  of  ferrite 
material,  TransTech  G-4259,  biased  with  a  1  kG  NDFE 
permanent  magnet.  The  thickness  of  the  disk  was  chosen 
to  be  0.344  inches,  producing  a  45  degree  Faraday  rotation 
of  the  microwave  beam  as  it  passes  through  the  disk.  Fused 
quartz  matching  layers  of  thickness  0.0425  inches  were  em¬ 
ployed  on  each  side  of  the  ferrite  disk  in  order  to  match  to 
free  space.  The  beam  was  then  injected  into  the  high  power 
overmoded  waveguide  system  via  an  identical  ellipsoid  and 
system  of  flat  mirrors. 


The  high  power  overmoded  waveguide  system  consists 
of  a  series  of  circular  waveguide  mode  converters  which 
changed  the  injected  beam  from  the  HEu  gaussian  mode 
to  the  TEo3  gyrotron  operating  mode.  Due  to  geomet¬ 
rical/mechanical  considerations  a  mode  converter  system 
similar  to  that  considered  by  Doane  [8]  was  employed.  The 
1.970  inch  diameter  output  waveguide  of  the  gyrotron  was 
nonlinearly  tapered  to  1.400  inches,  where  a  four-ripple 
TE03-TE02  mode  converter  was  inserted.  A  further  ta¬ 
per  to  a  1.00  inch  diameter  was  included  and  followed 
by  a  four-ripple  TEo2~TEoi  mode  converter.  The  conver¬ 
sion  efficiency  of  these  devices  each  exceeded  96%,  The 
guide  was  then  bent  through  an  angle  of  54.5®  over  an 
arc  length  of  20.0  inches  converting  the  TEqi  mode  to  the 
TMii  mode.  Since  the  bend  used  a  sinusoidal  curvature 
distribution  rather  than  a  constant  radius  of  curvature, 
an  efficiency  of  99.9%  was  achieved.  Finally,  a  corrugated 
TMii~IIEii  mode  converter  and  horn  assembly  was  used 
to  generate  a  gaussian  beam  with  the  correct  waist  size  to 
couple  into  the  quasi-optical  system  [9].  All  of  the  mode 


conversion  efficiencies  were  computed  and  verified  qualita¬ 
tively  by  far  field  pattern  measurements.  An  overmoded 
directional  coupler  was  also  included  in  the  waveguide  sys¬ 
tem,  providing  test  signals  for  a  crystal  detector  and  spec¬ 
trum  analyzer. 

The  gyrotron  output  followed  the  reverse  path  back  to 
the  Faraday  rotator  where  it  underwent  another  45  degree 
rotation.  Since  it  was  then  cross-polarized  to  the  grid,  the 
high  power  beam  was  transmitted  through  the  grid  and 
was  absorbed  by  a  free-space  load.  The  quasi-optic  circu¬ 
lator  system  including  the  mode  converters  was  measured 
to  have  an  insertion  loss  of  1  dB  and  an  isolation  of  greater 
than  25  dB.  The  isolation  was  measured  by  injecting  sig¬ 
nals  into  the  injection  signal  horn  and  into  the  TEoi~TMii 
mode  converter  using  a  mode  launcher  and  measuring  the 
transmission  coefficient  in  both  directions  using  an  HP8510 
ANA.  All  cross-polarization  effects  introduced  by  the  mir¬ 
rors,  horns,  and  waveguide  system  are  therefore  included  in 
the  measurement.  Waveguide  losses  reduced  the  injection 
signal  an  additional  1  dB. 


Fig.  4.  Block  diagram  of  the  experiment. 


A  small  portion  of  the  gyrotron  output  power  was  re¬ 
ceived  by  a  horn  antenna  located  near  the  free-space  load 
and  was  sent  to  the  input  port  of  a  balanced  mixer.  There 
it  was  combined  with  a  sample  of  the  injection  signal  from 
the  waveguide  coupler  located  at  the  TWTA  output.  The 
resulting  mixer  output  voltage  which  was  proportional  to 
the  instantaneous  phase  difference  between  the  output  and 
injection  signals,  was  fed  into  a  digitizing  oscilloscope.  The 
digitized  data  were  recorded  by  a  computer  and  stored  for 
later  analysis.  Crystal  detector  waveforms  and  spectrum 
analyzer  displays  were  also  obtained  from  a  branch  circuit 
in  the  gyrotron  output  waveguide, 

III.  Theory  Predictions:  Phase-locking  Equation 
AND  Demands  on  Gyrotron  Power  Supplies 


A.  Phase-locking  Equations 

In  order  to  obtain  the  phase-locking  condition  in  a  di¬ 
rect  manner,  we  start  from  the  gyrotron  operation  equa¬ 
tions  presented  in  reference  [10].  The  reduced  gyrotron 
equations,  which  describe  an  annular  electron  beam  inter¬ 
acting  with  a  circular  TEmn  mode  near  cutoff,  have  been 
described  in  a  number  of  papers  [11-15].  Here  we  simply 
write  down  the  pertinent  results.  For  operation  at  the  sth 
harmonic,  the  wave  equations  are 

^  =  ^^3r{F)+j  cosip  (2a) 

dtb  A 

—  =0-Igi(F)-—smiP  (2b) 

where  ^  +  igi  is  the  gain  functional  defined  by 

gif)  =  -i  r  dUiOiP’e'^). 

^  J  —  oo 


The  angle  brackets  denote  an  average  over  initial  phase  (j> 
uniformly  distributed  between  0  and  2t.  In  these  equa- 
tions, 

p=  — 

PlO 

is  the  perpendicular  momentum  normalized  to  its  initial 
value, 

(p=Wt-  S4>gp 

is  the  slowly  varying  gyrophase  {(i>gp  is  the  actual  particle 
gyrophase), 

2z 


L  ’ 


is  normalized  distance  where  L  is  the  length  of  the  cavity 
as  determined  by  the  Gaussian  axial  distribution, 

/(O  = 

and  /  is  the  normalized  current;  i.e.. 
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Fig.  5.  Quasi-optical  circulator  configuration. 


where  Tp  =  is  the  transit  time,  v^o  and  vxo  are 

the  axial  and  perpendicular  velocities,  respectively,  c  is  the 
speed  of  light,  is  the  initial  relativistic  cyclotron  fre* 
quency,  is  the  beam  current,  f q  is 

the  Alfven  current  in  mks  units  with  q  and  mo  the  electron 
charge  and  rest  mass,  respectively,  Q  is  the  loaded  quality 
factor  of  the  resonator,  kx  is  the  perpendicular  wave  num¬ 
ber  (^:x  ^  since  we  are  neeu*  cutoff),  Vg  is  the  guiding 
center  radius,  Vmn  is  the  nth  zero  of  the  derivative 
of  the  mth  order  Bessel  function,  and  m  is  the  azimuthal 
mode  number.  The  normalized  field  amplitude,  is 

given  in  terms  of  the  actual  field  as 


where  the  ‘hot’  resonant  frequency  of  the  cavity  is  uq  and 

Wo  Wp. 

Equation  (8)  together  with  Eqs.  (5)  and  (6)  indicates 
that  the  phase  traverses  180  degrees  as  the  signal  frequency 
swings  through  the  full  lock  bandwidth.  This  will  be  im¬ 
portant  in  demonstrating  that  phase-locking  was  accom¬ 
plished  during  the  experiments  by  injecting  an  external 
signal  that  is  swept  in  frequency. 

The  locking  bandwidth  defined  by  Eq.  (7)  coincides  with 
Adler’s  relation  [16].  Let  us  interpret  this  relation  in  terms 
of  the  input  and  output  power  flows  propagating  in  our 
system  in  the  same  waveguide  since  these  power  flows  can 


Bo  2vxo/c  si  I  2c  J  Jrn±>ikxrs) 
where  the  perpendicular  component  of  the  electric  field, 
Ex  =  Re  X  Vx«/m(i±r)exp  [i(m0  -  wi)]}  , 


be  measured  in  the  experiment.  For  this  purpose  let  us 
recall  that,  according  to  reference  [17],  the  flow  of  the  power 
into  the  cavity,  Fcav>  is  related  to  the  flow  of  the  incident 
wave  power,  P,n,  as 


^Qe/Qohm 


{l  +  Qe/QoUmV  +  A^^ 


(9) 


defines  Eq,  and  ^  is  the  slowly  varying  wave  phase  with 
respect  to  the  phase  of  the  injected  signal.  The  normalized 
frequency  mismatch,  /?,  is  given  by 


where  ujr  is  the  real  part  of  the  complex  cold  cavity  res¬ 
onant  frequency.  Its  imaginary  part,  which  is  responsible 
for  losses,  is  described  by  the  term  1/2  in  the  brackets  on 
the  left-hand  side  of  Eq.  (2a);  note  that  in  Eq.  (2),  r  is  the 
slowly  varying  time,  r  =  ut/Q,  Finally,  A  is  proportional 
to  the  field  injected  into  the  cavity  from  an  external  source. 
As  discussed  in  the  introduction,  this  paper  is  specifi¬ 
cally  addressed  to  phase-locking  of  a  harmonic  gyrotron 
with  external  injection.  In  this  respect,  the  dependance  of 
the  locking  bandwidth  A/x  on  the  input  and  output  power 
is  of  practical  interest.  This  quantity  can  be  estimated  by 
setting  d0/ dr  =  0  in  Eq.  (2b).  Then,  we  can  obtain  the 
following  formulae  for  the  maximum  and  minimum  detun¬ 
ing: 


/^max  —  ^9i  +  2^  Q  ^  sin  ^  =  1  (5) 

^min  =  l9i-^  =  Q  for  sin V-  =  -1  (6) 


Now,  if  we  assume  that  at  both  /?n,ax  and  the  imag¬ 
inary  part  of  the  gain  function,  gi,  and  the  normalized 
field  amplitudes  are  the  same,  then  we  can  obtain  the  fol¬ 
lowing  approximate  formulae  for  the  locking  bandwidth. 


^max  ~  ^min: 


^max  ^min  ^ 
CJQ  * 


(7) 


and,  for  the  dependence  of  phase  on  the  frequency: 


^(cj)  =  arcsin 


2(u;  —  cjo) 

' ) 

^max  “  ^min 


and  in  the  stationary  regime  this  power,  Pcav,  is  equal  to 
the  power  of  ohmic  losses  in  the  cavity, 


Tcav  —  P »hm  — 


Qohm 


A^V. 


(10) 


In  Eqs.  (9)  and  (10),  Qg  and  Qohm  aroj  respectively,  the  ex¬ 
ternal  and  ohmic  quality  factors  and  are  related  by  Q7*  = 
~  Qohm-  ^9-  (9)i  A  =  2Qg(oj  —  wo)/wo  and  in 
Eq.  (10)  V  is  the  resonator  volume.  At  the  same  time,  in 
the  presence  of  the  beam  when  the  microwave  power  ex¬ 
tracted  from  the  beam  is  much  larger  than  P-avi  the  output 
power  flow  is  equal  to 


e.u.  =  ~F’v.  (11) 

In  Eqs.  (10)  and  (11)  we  are  not  concerned  with  the  nor¬ 
malization  constant  for  A  and  F  since  they  are  the  same 
for  both,  and  we  are  interested  only  in  the  ratio  of  these 
values.  As  follows  from  Eqs.  (9)-(ll),  this  ratio  is  equal  to 


^  ^  2  I  Pm/ Pout  ' 

F  y  il+Qe/Qohn,r  +  A^- 


In  our  experiment,  the  condition  <  1  is  well  satisfied 
so  that  Eq.  (12)  reduces  to 

2  _ _ 

F  "  1  +  Qe/(?ohm  (13) 

Substituting  the  field  amplitude  ratio  given  by  Eq.  ( 13)  into 
Eq.  (7)  gives  the  following  equation  for  the  phase-locking 
bandwidth 


^max  ^ 


2 

<5(1  +  Qe/Qohm) 
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(14) 


For  the  case  of  a  critically  coupled  cavity  (i.e.,  Qg  =  Qohm) 
with  no  power  reflection  at  the  cavity  entrance,  Eq.  (14) 
becomes 
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(8) 


(15) 


Equation  (15)  is  another  common  form  of  Adler’s  equation; 
however,  it  is  not  applicable  to  the  present  case  since  the 
gyrotron  cavity  is  overcoupled  <  Qohm)-  Then,  one 
must  use  the  more  general  form  in  Eq.  (14),  which  can  be 
rewritten 

— ■j""'" = (16) 

We 

Equation  (16)  may  be  regarded  as  a  generalization  of  Adler’s 
equation  which  allows  for  power  reflection  at  the  entrance 
to  the  cavity. 

B,  Expected  Pulse  Waveforms 

In  the  experiment,  the  lock  range  can  be  determined  by 
sweeping  the  synthesizer  frequency  across  the  normal  op¬ 
erating  frequency  of  the  gyrotron  and  viewing  the  phase 
difference  between  the  gyrotron  output  and  the  injected 
signal  over  a  single  pulse  as  well  as  on  a  pulse-to-pulse  ba¬ 
sis. 


•  MIXER  VOLTAGE  (PHASE  DIFFERENCE)  (AVERAGED  OVER  PULSE  ) 

«  STANDARD  DEVIATTON  OF  THE  MIXER  VOLTAGE  OVER  PULSE 
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Fig.  6.  Expected  phase  vs,  injected  frequency  over  three 
regions:  1)  high  frequency  beat  note,  2)  uncorrelated  phase 
difference,  3)  phase-locked  for  the  specific  case  where  the 
180®  range  of  phase  difference  is  from  0®  to  180®. 


The  average  value  of  the  phase  difference  between  the 
reference  and  test  signals  over  a  given  pulse  is  estimated 
by  computing  the  average  of  a  set  of  samples  taken  from  the 
center  portion  of  the  difference  pulse.  When  the  difference 
frequency  is  large,  many  cycles  appear  across  the  pulse  and 
the  average  voltage  is  zero.  When  the  two  input  frequen¬ 
cies  to  the  mixer  are  close  enough  so  that  no  significant 
portion  of  a  cycle  appears  across  the  pulse,  a  finite  average 
mixer  voltage  is  found  for  that  particular  pulse.  If  the  in¬ 
jection  signal  is  so  close  to  the  gyrotron  frequency  that  the 
difference  is  not  resolvable  over  the  pulse  length,  but  the 
gyrotron  is  not  phase-locked  with  the  injection  signal  we 
would  expect  uncorrelated  phase  (average  pulse  voltage) 
from  pulse  to  pulse.  When  the  gyrotron  is  injection  locked 
we  expect  the  phase  to  ramp  through  a  180  degree  range  as 
the  locking  bandwidth  is  traversed  (see  the  next  section). 
An  example  of  average  voltage  results  over  three  regions 
(high  frequency  beat  note,  uncorrelated  phase  difference, 
and  phase-locked)  are  shown  in  Fig.  6.  In  addition  to  the 


average  phase  of  the  pulse,  the  standard  deviation  of  the 
mixer  voltage  (pulse  phase)  is  also  included  on  the  plot. 
This  statistic  is  also  an  aid  in  determining  the  lock  range. 
When  many  cycles  appear  across  the  pulse  the  standard 
deviation  is  high;  it  is  low  under  locked  conditions  regard¬ 
less  of  the  phase  of  the  pulse. 

(7.  Demands  on  Gyrotron  Beam  Voltage 

In  real  systems,  there  are  always  voltage  fluctuations 
which  change  the  hot  cavity  resonant  frequency,  cjq*  For 
phase-locking  to  occur,  frequency  changes  caused  by  the 
voltage  fluctuations,  which  we  denote  must  be  smaller 
than  the  locking  bandwidth,  Combining  this  condi¬ 

tion  with  Eq.  (16)  yields 


< 


2ujo 


(17a) 


Writing  Aw^  =  AVidcjo/dH,  where  A  Vi  is  the  voltage 
fluctuation,  we  have 


^Vb  j  Pm  2uJq _ 

Vi  VPout  Vige(da;o/dVi)' 


(17b) 


In  this  equation,  A  Vi/ Vi  is  to  be  interpreted  as  the  al¬ 
lowed  relative  voltage  fluctuation  in  the  modulator.  Thus, 
if  Vic/cJo/dVi  is  too  large,  it  will  be  impossible  to  construct 
a  voltage  modulator  with  small  enough  fluctuations  to  sat¬ 
isfy  Eq.  (17).  We  reproduce  here  the  calculated  results 
which  appeared  in  reference  [10]  and  which  are  shown  in 
Fig.  7.  Allowed  maximum  values  of  AVi/Vi  are  plotted 
versus  voltage  for  various  values  of  the  detuning  parame¬ 
ter  8  which  is  given  by 
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where  Qco  is  the  initial  relativistic  cyclotron  frequency. 

In  these  figures,  the  assumption  was  made  that  Pm/ Pout  = 
10''^  and  that  Qc  =  1000.  The  value  of  and  the  range 
of  Vb  do  not  match  those  in  the  present  experiment,  but 
the  trend  of  the  predictions  are  believed  to  be  relevant. 
Gyrotron  efficiency  is  also  included  on  the  plots.  In  the 
region  where  the  efficiency  is  below  approximately  20%, 
the  requirement  on  the  modulator  voltage  fluctuations  are 
quite  stringent:  on  the  order  of  0.1%.  On  the  other  hand, 
for  gyrotron  efficiency  ^  30%,  the  allowed  voltage  fluctua¬ 
tions  are  typically  on  the  order  of  1%. 


Beam  VolUge  (kV) 
c)  (5=0.35 


'  Allowed  voltage  fluctuation  for  phase-locking. 


IV.  Experimental  Results 


A,  Measurement  of  Locking  Bandwidth 

In  all,  sixteen  separate  experiments  were  performed.  Two 
were  control  experiments  with  the  injection  horn  blocked, 
eight  were  performed  at  full  injection  power  (60  Watts)! 
and  two  were  performed  at  each  of  three  reduced  power 
levels,  (30,  15,  and  6  Watts).  For  each  experiment  250 
pulses  were  captured  near  the  center  of  the  lock  range  for 
plater  analysis.  The  synthesizer  was  swept  across  a  10  MHz 
range  in  60  seconds,  and  the  pulse  repetition  frequency 
was  43  Hz,  corresponding  to  a  frequency  step  of  3.876  kHz 
per  pulse.  The  pulse  waveforms  were  sampled  at  a  rate  of 
10  samples  per  microsecond,  and  the  effective  pulse  width 
was  4  microseconds.  In  some  instances  the  250  pulses  did 
not  straddle  the  complete  locking  bandwidth  and  no  useful 
results  were  obtained  from  those  experiments. 

Figure  8  plots  the  average  mixer  voltage  (proportional  to 
the  sine  of  the  phase),  using  a  solid  line  to  connect  each  of 
the  250  points  in  a  control  experiment  where  the  injection 
horn  was  blocked.  As  expected,  near  the  beginning  and 
near  the  end  of  the  data  the  difference  frequency  is  large 
enough  to  produce  at  least  one  complete  cycle  across  the 
pulse  and  zero  average  voltage  is  obtained.  Between  pulse 
50  and  100  the  difference  frequency  is  small,  and  only  the 
random  pulse-to-pulse  phase  jitter  of  the  gyrotron  is  seen. 
No  pulse-to-pulse  correlation  is  observed  and  the  average 
voltage  fluctuates  wildly  across  the  entire  allowable  mLxer 
range  of  0.35  to  -0.25  volts.  Note  the  fact  that  these  two 
voltage  extremes  are  not  symmetric,  indicating  a  non-ideal 
mixer  with  some  finite  offset  voltage.  This  plot  is  exactly 
what  one  would  expect  for  the  control  experiment  with 
no  injection  power.  The  standard  deviation  of  the  voltage 
samples  on  each  pulse  is  plotted  as  a  dashed  curve  in  Fig.  8. 
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Fig.  8.  Mean  and  standard  deviation  of  samples  (no  lock¬ 
ing)- 

Figure  9  plots  the  same  information  for  the  case  when 
the  termination  is  removed  and  60  Watts  are  injected  into 
the  gyrotron.  Now  a  constant  phase  progression  covering 
the  full  mixer  range  (180  degrees)  is  seen  from  pulse  num¬ 
ber  50  to  164.  This  114  pulse  interval  then  represents  the 
locking  bandwidth,  which  in  this  case  corresponds  to  440 
kHz.  It  can  also  be  observed  that  at  the  high  end  of  the 
lock  range,  near  pulse  number  160,  there  is  a  region  where 
intermittent  locking  occurs.  Shot-to-shot  modulator  pulse 
height  fluctuations  cause  shifts  in  the  gyrotron  output  fre¬ 
quency  sufficient  to  randomly  unlock  the  system.  This  can 
be  attributed  to  the  limited  pulse-to-pulse  stability  as  well 
as  the  long  term  drift  of  the  gyrotron  frequency.  Including 
the  intermittent  region,  locking  occurs  over  a  bandwidth 
of  510  kHz. 
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Fig.  9.  Mean  and  standard  deviation  of  samples  (locked 
114  pulses  -  440  kHz). 


Examination  of  the  individual  pulses  in  the  locking  range 
shows  that  the  behavior  predicted  in  Fig.  6  is  observed. 


Figure  10  plots  the  mixer  voltage  for  four  representative 
pulses  captured  during  the  experiment  discussed  above.  As 
can  be  seen  from  Fig.  10,  pulse  number  3  is  unlocked  and 
more  than  a  complete  cycle  of  beat  note  is  visible  across  the 
pulse.  Pulses  numbered  75,  120,  and  140  represent  locked 
pulses.  They  are  relatively  flat,  with  an  average  voltage 
that  is  proportional  to  the  phase  difference  between  the 
injection  and  output  signals.  This  voltage  varies  from  a 
maximum  (pulse  75)  through  zero  (near  pulse  120)  to  a 
minimum  (pulse  140)  as  the  lock  range  is  traversed. 

In  all,  twelve  of  the  experiments  produced  useful  results. 
These  data  are  tabulated  below  in  Table  3  along  with  the¬ 
oretical  predictions  from  Eq.  (16)  with  gyrotron  output 
power  Pout  =  100  kW,  and  the  gyrotron  operating  fre¬ 
quency  /o  =  34.5  GHz,  and  taking  =  2500.  There  is 
rough  agreement  between  the  experimental  values  of  lock¬ 
ing  bandwidth  and  the  theoretical  predictions.  At  the  high¬ 
est  values  of  the  injection  power,  the  measured  bandwidth 
is  somewhat  smaller  than  the  theoretical  value;  this  may 
be  due  to  the  ripple  in  the  gyrotron  modulator  voltage.  We 
also  note  that  as  injection  power,  P,*,  is  reduced  the  locking 
bandwidth  appears  to  shrink  less  rapidly  than  \fPi\  how¬ 
ever,  further  study  is  required  before  a  definitive  statement 
can  be  made. 
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Fig.  10.  Mixer  voltage  for  four  representative  pulses  cap¬ 
tured  during  the  experiment. 


Table  3.  Results  of  phase-locking  measurement. 


P.n 

(Watts) 

Avg.  # 
Pulses 
Locked 

Bandwidth 

(kllz) 

Theoretical 
Bandwidth 
(kHz)  using 
Eq.  (17) 

Gain 

(dB) 

0 

0 

0 

0 

0 

6 

68 

230-310 

220 

42 

15 

79 

310-330 

348 

38 

30 

90 

340-390 

490 

35 

60 

106 

440-510 

693 

32 

B.  Effect  of  Beam  Voltage  on  Gyrotron  Operation 
The  phase-locking  experiment  was  performed  with  a  gy¬ 
rotron  operating  efficiency  of  30%  and  with  a  beam  voltage 
ripple  of  0.3%.  These  conditions  are  compatible  with  theo¬ 
retical  estimates  for  the  realization  of  phase-locking  as  sug¬ 
gested  by  Fig.  7.  In  fact,  the  pulse  width  of  the  beam  volt¬ 
age  was  7  microseconds,  but  the  duration  of  the  microwave 
signal  over  which  phase-locking  could  be  realized  was  4  mi¬ 
croseconds.  This  4  microsecond  interval  corresponds  to  the 
flat  portion  of  the  beam  voltage  where  AV^/Vj  <  0.3%. 

V.  Discussion  of  Results 

The  present  experiment  demonstrated  phase-locking  of  a 
second-harmonic,  100  kW,  millimeter-wave  gyrotron  oscil¬ 
lator.  Locking  bandwidth  was  measured  by  means  of  an  in¬ 
stantaneously  recording  and  computer  processing  system. 
The  results  demonstrate  for  the  first  time  phase-locking  of 
a  harmonic  gyrotron  oscillator  (s  >  1);  the  phase-locking 
was  accomplished  by  separating  the  paths  of  the  injection 
signal  and  the  output  signal  in  a  low  loss,  high  power  com¬ 
patible,  quasi-optical  circulator.  The  insights  gained  dur¬ 
ing  this  study  may  be  summarized  as  follows: 

(a)  A  second-harmonic  gyrotron  oscillator  incorporating 
a  special  complex  cavity  can  operate  with  sufficient 
stability  and  efficiency  that  it  can  be  phase-locked  by 
injecting  an  external  signal  as  small  as  42  dB  below 
the  gyrotron  output  power. 

(b)  A  quasi-optical  circulator  can  be  of  practical  value  in 
phase-locking  high  power  millimeter-wave  oscillators 
with  signal  injection  through  the  output  window. 

(c)  The  requirement  for  gyrotron  power  supply  stability 
to  achieve  phase-locking  is  quite  demanding  but  can 
be  realized  if  gyrotron  efficiency  is  high  (i.e.,  AVj/I4  ~ 
0.3%  is  required  if  efficiency  is  ^  30%). 

(d)  Measured  values  of  locking  bandwidth  roughly  agree 
with  the  predictions  of  a  generalized  Adler’s  equation. 
The  measured  values  of  locking  bandwidth  for  a  given 
value  of  the  ratio,  P\n/Po\xty  were  much  smaller  in  the 
present  study  than  in  a  previous  study  of  a  gyrotron 
operating  at  the  fundamental  of  the  electron  cyclotron 
frequency  [5];  the  difference  is  expected,  and  is  at¬ 
tributable  to  the  longer  cavity  length  and  correspond¬ 
ingly  larger  value  of  in  the  gyrotron  of  the  present 
study. 
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Abstract 

We  present  the  theoretical  design  of  a  second  harmonic  small-orbit  gyrotron  ampli¬ 
fier  which  utilizes  the  interaction  between  a  35  kV,  4  A  beam  and  a  TEqh  cavity  to 
produce  over  70  kW  of  amplified  power  at  9.9  GHz  in  a  1.83  kG  magnetic  field.  One  of 
the  novel  features  of  this  device  is  that  the  electron  gun  produces  an  axially-streaming 
annular  beam  which  is  velocity  modulated  by  a  short  TMono  input  cavity.  Perpendicular 
energy  is  imparted  to  the  beam  via  a  non-adiabatic  magnetic  transition  at  the  end  of  a  13 
cm  drift  region.  An  efficiency  of  53%  is  predicted  with  a  large  signal  gain  near  20  dB 
by  a  single  particle  code  which  takes  into  account  non-ideal  effects  associated  with  fi¬ 
nite  beam  thickness  and  finite  magnetic  field  transition  widths.  A  description  of  the 
operating  principles  and  details  of  the  modeling  techniques  for  this  device  are  given  be¬ 
fore  the  tube  performance  is  characterized. 


I.  INTRODUCTION 


In  recent  years,  numerous  theoretical  and  experimental  investigations  of  gyrotrons  operating 
near  the  fundamental  cyclotron  frequency  have  proven  them  to  be  reliable,  efficient,  high  power 
sources  of  microwave  and  millimeter  wave  radiation.''"^  These  studies  have  been  performed  over 
a  wide  range  of  parameters  and  have  resulted  in  a  wealth  of  information  about  the  electron  gun 
requirements  and  characteristics  and  about  the  stability,  efficiency,  and  mode  competition  prop¬ 
erties  of  microwave  circuits  for  gyrotrons.^’^  Potential  applications  of  these  sources  include  RF 
drivers  for  accelerators,  plasma  heating  in  magnetic  fusion  research,  millimeter  wave  and  deep 
space  radar,  materials  processing,  and  nonlinear  spectroscopy  of  semiconductors  and  biological 
materials. 

Unfortunately,  conventional  gyrotrons  operating  near  the  fundamental  cyclotron  frequency 
(first  harmonic  gyrotrons)  must  be  immersed  in  a  magnetic  field  whose  strength  is  linearly  pro¬ 
portional  to  the  output  frequency.  Thus,  they  are  not  viable  candidates  for  some  high  frequency 
applications  that  require  compact  and  lightweight  tubes.  Two  potential  ways  to  decrease  the  ap¬ 
plied  magnetic  field  and  the  accompanying  material  bulk  and  power  consumption  are  (1)  to  use 
Doppler  upshifting  via  operation  at  a  high  axial  wavenumber  and  (2)  to  operate  at  a  harmonic  of 
the  cyclotron  frequency.  The  cyclotron  auto-resonance  maser  (CARM)  is  an  example  of  a  device 
that  takes  the  first  approach.  Many  CARMs  have  been  built  and  tested  but  none  have  come  close 
to  the  efficiencies  achieved  in  first  harmonic  gyrotrons.*^''^  Mode  competition  from  gyrotron  in¬ 
stabilities  and  sensitivity  to  velocity  spread  have  typically  been  the  limiting  factors.  Low 
harmonic  operation  (typically  the  second  or  third  harmonic)  has  met  with  some  success  in  the 
conventional  small-orbit  gyrotron  configuration'^’'"*  and  high-harmonic  devices  have  been  fairly 
successful  in  the  large-orbit  configuration. Still,  the  efficiencies  of  second  harmonic  devices 
usually  lag  behind  their  first  harmonic  counterparts  by  several  percent  and  efficiencies  tend  to 
decrease  fairly  rapidly  with  increasing  harmonic  number  after  that  point. 
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We  have  recently  proposed^'*  a  method  of  pre-bunching  the  electron  beam  to  enhance  the  ef¬ 
ficiency  in  high  harmonic,  large-orbit  devices  which  is  based  on  the  cusp-injection  scheme.^^ 
Before  the  annular  beam  passes  through  the  balanced,  non-adiabatic  magnetic  field  reversal,  it 
first  encounters  a  circularly  polarized  TM„,„o  mode  input  cavity  and  a  subsequent  drift  region. 
The  resultant  ballistic  bunching  (potentially)  sets  up  a  beam  that  can  efficiently  interact  at  the 
ffjth  harmonic  with  a  simple  right-circular  output  cavity.  For  example,  the  design  study  in  Ref  34 
demonstrated  that  an  efficiency  of  40%  was  achievable  in  X-Band  at  the  fourth  harmonic  via  an 
interaction  between  a  100  kV,  25 A  electron  beam  and  a  TE41 1  cavity. 

While  this  initial  result  is  quite  promising,  there  are  some  practical  limits  to  the  applicable 
range  of  parameters  for  this  device.  For  example,  the  beam  radius  of  the  large-orbit  configuration 
is  equal  to  the  Larmor  radius,  which  is  proportional  to  the  perpendicular  velocity  for  a  given  cy¬ 
clotron  frequency.  Thus,  the  beam  radius  in  tow  voltage  systems  decreases  fairly  substantially 
with  decreasing  beam  energy.  However,  the  relative  radius  of  the  maximum  electric  field  in  a 
cavity  increases  with  increasing  harmonic  number.  In  fact,  the  electric  field  at  radii  much  less 
than  the  maximum  field  radius  varies  by  a  factor  that  is  proportional  to  the  power  of  the  har¬ 
monic  number.  In  such  systems,  the  achievement  of  the  necessary  field  strength  for  efficient 
interaction  often  leads  to  unrealistic  requirements  on  output  cavity  quality  factors. 

There  are  at  least  two  possible  ways  to  overcome  these  difficulties  and  design  efficient,  low 
voltage,  harmonic  devices.  One  way  is  to  use  a  vane-resonator  to  couple  various  azimuthal 
modes  together.  The  resultant  cavity  modes  can  support  the  necessary  azimuthal  field  variation 
with  a  much  lower  cutoff  frequency.  Thus,  the  wall  radius  can  be  reduced  considerably  as  com¬ 
pared  to  the  smooth  wall  case  and  the  interaction  impedance  can  be  significantly  increased.  There 
are  some  potential  drawbacks  to  this  approach.  For  example,  the  mode  composition  of  the  output 
signal  in  an  axially-coupled  vane  resonator  cavity  can  be  significantly  more  complicated  than  for 
the  smooth  wall  case.  Still,  this  approach  has  been  used  with  reasonable  success  in  many  con- 
ventional  large-orbit  gyrotrons.  ^  The  design  of  an  axially-bunched,  low  voltage  large-orbit 
gyrotron  with  a  vane  resonator  will  be  the  subject  of  a  later  study. 
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The  approach  we  will  investigate  in  this  work  is  to  pre-bunch  a  small-orbit  beam.  Axial 
modulation  is  still  a  viable  concept  because  a  small  orbit  beam  can  be  generated  by  a  non- 
adiabatic  magnetic  transition  that  does  not  reverse  the  direction  of  the  field.  The  guiding  center 
radius  can  be  adjusted  to  place  the  beam  at  a  location  where  the  electric  field  is  relatively  strong. 
This  non-zero  guiding  center  radius  expands  the  number  of  potential  modes  that  the  beam  can 
interact  with.  While  this  may  enable  us  to  select  a  desirable  output  mode  it  also  can  potentially 
lead  to  additional  mode  competition  problems.  A  properly  pre-bunched  beam,  however,  should 
help  to  reduce  mode  competition  by  preferentially  selecting  the  desired  operating  mode. 

We  will  illustrate  this  approach  via  the  presentation  of  a  second  harmonic  small-orbit  system 
design:  In  this  design,  a  35  kV,  4  A  beam  is  pre-bunched  by  a  TM020  cavity  and  interacts  with  a 
TEoii  cavity  to  theoretically  produce  over  70  kW  of  power  with  an  efficiency  and  large-signal 
gain  of  approximately  53%  and  20  dB,  respectively.  In  Section  II,  we  discuss  the  basic  amplifier 
configuration  and  describe  the  computer  codes  used  in  the  analysis.  In  Section  III,  we  detail  the 
results  of  the  design  example  and  present  tube  stability  and  large-signal  characteristics.  We  also 
characterize  the  sensitivity  of  the  device  to  parameter  variations.  The  results  of  this  study  are 
summarized  in  the  final  section. 
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II.  THEORETICAL  MODEL 


A  simplified  schematic  of  an  axially-modulated,  cusp-injected,  small-orbit  gyrotron  ampli¬ 
fier  is  given  in  Fig.  la.  A  linearly-streaming,  thin,  mono-energetic,  annular  beam  is  immersed  in 
a  uniform  axial  magnetic  field  and  injected  into  a  short  buncher  cavity.  The  axial  electric  field  of 
the  input  cavity  (Fig.  Ic)  modulates  the  energy  of  the  electron  beam  which  results  in  ballistic 
bunching  as  the  electrons  drift  through  the  region  immediately  following  the  cavity.  When  the 
electrons  encounter  the  non-adiabatic  magnetic  field  transition  (shown  in  Fig.  lb),  the  radial 
magnetic  field  produces  a  Lorentz  force  which  converts  part  of  each  electron’s  energy  into  per¬ 
pendicular  motion.  Because  this  conversion  occurs  in  a  short  distance,  the  electron  bunches  that 
form  begin  to  rotate  with  the  particles.  These  bunches  then  interact  with  the  electric  field  in  the 
output  cavity  (Fig.  Ic)  and,  if  properly  formed,  can  give  up  most  of  their  energy  to  the  electro¬ 
magnetic  wave. 

In  order  for  the  bunches  to  interact  effectively  with  the  output  cavity  fields,  the  resonant  fre¬ 
quencies  of  the  input  and  output  cavities  must  be  nearly  integer  multiples  of  the  cyclotron 
frequency  (in  the  output  cavity).  Furthermore,  if  Sj  and  S2  are  the  multiples  (harmonics)  of  the 
cyclotron  frequency  in  the  input  and  output  cavities,  respectively,  and  mj  and  m2  are  the  azi¬ 
muthal  mode  indices  of  the  same  two  cavities,  respectively,  then  the  following  condition  must  be 
satisfied: 

mi/si=m2/s2.  (1) 

The  designs  in  this  work  utilize  operating  modes  with  m,  =  m2  =  0,  so  that  any  combination  of  s, 
and  S2  are  acceptable.  We  will  always  focus  on  the  S2  =2  case.  Most  of  the  analysis  will  be  per¬ 
formed  with  Si  =2  as  well,  but  we  will  demonstrate  towards  the  end  that  an  Si=l  design  is  nearly 
as  effective  as  the  earlier  design  and  has  several  nice  features  that  warrant  consideration. 

Because  the  input  cavity  operates  in  the  TMono  mode,  it  can  be  excited  either  by  a  single  ap¬ 
erture  on  the  radial  wall  fed  by  a  standard  waveguide  or  by  a  single  loop  fed  by  a  simple  coaxial 
line.  The  required  drift  tube  length  d  depends  on  the  magnitude  of  the  energy  spread  {AE/E^  im¬ 
parted  to  the  beam  and  has  been  derived  elsewhere  for  a  low  voltage  system  as:^"* 
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d=^,%l{2AE/E,)  (2) 

where  is  the  streaming  velocity  (normalized  to  the  speed  of  light)  and  X  is  the  RF  wavelength 
(of  the  input  cavity  signal). 

The  magnetic  transition  at  the  end  of  the  drift  region  is  shortened  by  a  thin  iron  pole  piece. 
Long  solenoids  on  either  side  of  the  pole  piece  generate  approximately  constant  axial  magnetic 
fields  in  each  region  (B^  on  the  input  side  and  Bg  on  the  output  side).  The  required  mismatch 
between  the  magnetic  field  strengths  depends  on  the  beam  radius  and  the  desired  velocity  ratio 
If  the  original  beam  radius  is  r^  and  the  Larmor  radius  after  the  magnetic  transition  is 

or 

rL,  then  the  ratio  for  an  ideal  (zero  length)  transition  is: 

Be/Bo=l-2rL/r„.  (3) 

The  beam’s  energy  spread  results  in  velocity  spreads  after  the  magnetic  transition  which  can 
be  approximated  by:^"* 


Av 


z  _ 


1  +  a^  +a^(/l-l) 


AE 


^1.  2  E^ 


(4) 


when  the  effects  of  the  input  cavity’s  magnetic  fields  are  neglected.  In  Eq.  (4),  Yq  is  the  relativis¬ 
tic  mass  factor.  For  low  voltage  beams,  the  induced  perpendicular  velocity  spread  is  typically 
small  but  the  parallel  spread  can  become  quite  large.  Setting  an  upper  limit  on  this  value  often 
dictates  the  drift  tube  length  in  small-orbit  designs  via  Eq.  (2). 

There  are  two  other  contributions  to  the  net  velocity  spreads.  The  first  is  from  the  electron 
gun  optics  and  hence  is  design  specific,  can  typically  be  kept  fairly  small,  and  won’t  be  discussed 
further.  The  other  contribution  comes  from  the  canonical  angular  momentum  of  finite  thickness 
beams.  This  spread  can  be  reduced  in  small-orbit  systems  by  compressing  the  beam  width.^^ 

The  output  cavity  for  the  pre-bunched  small  orbit  gyrotron  is  a  simple  axially-coupled  right- 
circular  waveguide  that  is  designed  to  operate  at  the  proper  frequency  in  the  TEqh  mode.  The 
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cavity  has  a  long  main  section  and  a  shorter  output  lip  that  has  a  radius  smaller  than  that  required 
for  the  operating  mode  to  propagate.  The  resonant  frequency  and  quality  factor  (Q)  are  calculated 
with  a  scattering-matrix  code.^’’^^  In  general,  the  required  cavity  Q  is  far  below  the  resistive  Q  of 
a  copper  cavity  but  sufficiently  large  so  that  the  interaction  between  the  beam  and  the  micro- 
waves  occurs  predominantly  in  the  main  section  of  the  cavity.  For  this  reason,  the  sinusoidal 
axial  field  dependence  of  a  closed  cavity  is  often  a  good  approximation  for  the  open-ended  cav¬ 
ity. 

Four  main  computer  codes  are  utilized  in  the  design  process.  The  scattering  matrix  code  is 
used  to  design  the  output  cavity  and  a  start-oscillation  code^^  is  used  to  check  the  stability  of  the 
output  cavity  to  spurious  modes  in  the  absence  of  beam  pre-bunching  by  the  input  cavity.  This 
code  assumes  that  the  magnetic  field  is  constant  in  the  output  cavity. 

The  other  two  codes  are  single  particle  codes  that  numerically  integrate  the  particle  motion 
in  the  device.  The  first  single  particle  code  simulates  the  beam’s  motion  through  the  input  cavity, 
the  drift  region,  and  the  non-adiabatic  magnetic  transition.  Given  the  beam  parameters,  electrons 
are  initially  “launched”  over  a  range  of  axial  positions,  times,  and  angles  that  represent  the  en¬ 
tirety  of  phase  space  (after  taking  into  account  any  symmetries  in  the  system).  The  initial 
velocities  of  the  particles  are  assumed  to  be  purely  axial,  but  electrons  can  be  started  from  several 
radial  locations  to  emulate  finite  thickness  beams.  The  EM  fields  are  assumed  to  be  confined  to 
the  input  cavity  (leakage  fields  are  neglected).  The  only  field  in  the  drift  region  is  the  axial  mag¬ 
netic  field,  which  is  assumed  to  remain  constant.  The  magnetic  transition  region  is  simulated 
with  a  piece-wise  linear  axial  field  profile.  The  Runga-Kutta  method  is  used  to  start  the  process 
and  a  predictor-corrector  scheme  is  used  for  the  bulk  of  the  numerical  integration.  Numerical 
convergence  is  checked  by  varying  integration  step  size  and  the  number  of  particles.  Conserva¬ 
tion  of  energy  and  canonical  angular  momentum  are  also  checked  in  regions  where  there  is  only 
a  static  magnetic  field.  This  code  produces  the  phase  space  distributions  of  the  electrons  at  the 
end  of  the  transition  region  and  calculates  various  quantities  that  give  the  average  beam  statistics 
and  indicate  the  effectiveness  of  the  bunching  process.  The  resistive  Q  of  a  copper  input  cavity 
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and  the  required  drive  power  are  also  reported.  This  code  is  typically  iterated,  while  adjusting  the 
drift  region  length  and  the  strength  of  the  input  field,  until  satisfactory  bunching  performance  is 
achieved. 

The  second  code  takes  the  results  from  the  first  single  particle  code  and  uses  the  same  nu¬ 
merical  techniques  to  calculate  the  extraction  efficiency  of  the  EM  wave  in  the  output  cavity 
under  steady-state  conditions.  The  code  can  model  either  a  closed  circular  cavity  or  it  can  accept 
the  output  from  the  scattering  matrix  code  to  model  open  resonator  configurations.  The  magnetic 
field  in  the  output  cavity  can  be  tapered  linearly  to  optimize  the  interaction.  This  code  is  typically 
iterated,  while  varying  the  cavity  and  magnetic  field  parameters,  until  an  optimal  design  is 
achieved.  During  each  efficiency  run,  the  amplitude  and  phase  of  the  EM  wave  are  adjusted  until 
maximum  efficiency  is  calculated.  In  addition  to  the  efficiency,  the  code  reports  the  required  dif¬ 
fractive  Q  and  the  resistive  Q  for  a  closed,  right-circular  copper  cavity.  By  decreasing  the  field 
amplitude  in  the  output  cavity  toward  zero,  this  code  can  also  be  used  to  calculate  the  start  cur¬ 
rents  for  the  various  modes  in  tapered  magnetic  fields,  whether  the  beam  is  pre-bunched  or 
unbunched. 

We  have  also  run  a  separate  large  signal  code  to  independently  verify  the  efficiency  calcula- 
tions  in  the  output  cavity.  This  code,  originally  designed  for  relativistic  gyroklystrons  ,  is 
modified  to  take  into  account  the  pre-bunched  beam  at  the  entrance  of  the  output  cavity  from  the 
first  single  particle  code.  The  code  also  assumes  the  cold  cavity  field  profiles  and  Q  value  from 
the  scattering  matrix  code  which  includes  linear  mode  conversion  at  the  radial  steps  in  the  output 
cavity.  In  order  to  optimize  the  device  efficiency,  a  two-dimensional  search,  in  terms  of  overall 
mode  amplitude  and  phase,  is  carried  out.  To  ensure  self-consistency  in  the  calculations  the 
quality  factor  required  to  maintain  the  steady-state  field  amplitude  is  matched  with  the  Q  value 
from  the  scattering  matrix  code. 

Details  of  the  results  from  these  simulation  codes  for  the  specific  designs  are  given  in  the 
following  section. 
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III.  DESIGN  SIMULATION 


The  parameters  for  the  second  harmonic  design  are  given  in  Table  1.1.  The  system  is  de¬ 
signed  to  produce  microwaves  at  9.9  GHz  via  the  interaction  between  a  35  kV,  4  A  beam  and  the 
microwave  circuit.  The  initial  beam  radius  of  1.5  cm  is  selected  to  optimize  the  interaction  with 
the  EM  wave  in  the  output  cavity  while  avoiding  any  beam  interception.  This  beam  can  be  gen¬ 
erated  with  a  standard  electrode  configuration.^'*  If  a  magnetic  compression  of  four  were  to  be 
utilized,  for  example,  the  cathode  thickness  would  be  2  mm,  the  average  emitter  radius  would  be 
3  cm,  and  the  required  current  density  would  be  quite  modest  at  a  little  over  1  A/cm  .  The  table 
indicates  the  pre-  and  post-compression  beam  dimensions.  The  DC  space  charge  depression  in 
the  output  cavity  is  less  than  3  %  and  is  considerably  smaller  before  the  magnetic  transition.  The 
radial  thickness  corresponds  to  an  ideal  post-cusp  velocity  spread  of  Av_l=  1.9  %  (if  beam  width 
compression  is  not  utilized  to  reduce  this  number).  The  average  velocity  ratio  is  somewhat  larger 
than  a  typical  small-orbit  system  but  is  below  what  is  often  quoted  for  cusp-injected  large-orbit 
systems. 

The  drive  cavity  is  a  simple  right-circular  TM020  cavity  with  the  dimensions  indicated  in  the 
table.  The  cavity  width  is  selected  to  give  a  beam  coupling  factor'*®  of  about  0.93.  The  cavity  ra¬ 
dius  is  adjusted  to  achieve  the  desired  operating  frequency.  Note  that  the  quality  factor  quoted 
includes  only  the  losses  in  the  copper.  The  actual  quality  factor,  after  taking  into  account  the  in¬ 
put  power  coupling,  should  be  about  one-half  of  this  value.  All  drive  powers  quoted  in  this  paper 
assume  a  net  input  Q  of  1400.  Any  beam  loading  is  neglected.  The  beam  is  located  between  the 
first  null  and  the  second  peak  of  the  axial  electric  field.  The  required  input  power  could  be  de¬ 
creased  and  the  resultant  gain  increased  if  the  beam  could  be  re-positioned  on  the  maximum.  The 
easiest  way  to  accomplish  this  would  undoubtedly  be  to  use  a  coaxial  input  cavity.  Other  possi¬ 
bilities  include  the  introduction  of  dielectrics  into  the  input  cavity  or  moving  the  input  cavity  into 
the  magnetic  compression  region  where  the  beam  radius  is  larger. 

The  peak  axial  magnetic  field  in  the  output  cavity  is  determined  from  the  optimization  pro¬ 
cedure  and  is  such  that  twice  the  cyclotron  frequency  is  about  3%  lower  than  the  drive  frequency. 
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The  magnetic  field  in  the  input  region  is  about  60  Gauss  less  than  predicted  by  Eq.  (3)  due  to  the 
finite  cusp  width.  The  distance  between  the  input  cavity  and  cusp  transition  is  optimized  with  the 
numerical  codes,  but  the  initial  length  is  chosen  from  Eqs.  (2)  and  (4)  so  that  maximal  bunching 
could  be  achieved  with  an  energy  spread  and  an  induced  axial  velocity  spread  approximately  un¬ 
der  3%  and  10%,  respectively.  The  transition  width  should  readily  achievable  with  the  aid  of  one 
or  more  iron  pole  pieces.  The  optimal  magnetic  field  profile  is  shown  in  Fig.  lb. 

The  length  of  the  output  cavity  is  varied  to  maximize  efficiency  while  keeping  the  required 
Q  less  than  5-7%  of  the  resistive  Q  of  a  copper  cavity.  This  restriction  means  that  over  93%  of 
the  microwave  power  can  be  extracted.  The  cavity  radius  and  lip  dimensions  are  adjusted  to  pro¬ 
duce  the  required  resonant  frequency  and  Q  while  avoiding  beam  interception  (by  the  particles 
that  are  accelerated  during  the  interaction).  The  main  cavity  length  is  such  that  particles  undergo 
approximately  eleven  revolutions  inside  this  region.  The  azimuthal  electric  field  at  the  beam  ra¬ 
dius  in  the  output  cavity  is  shown  in  Fig.  Ic.  The  exponential  decay  in  the  drift  region  and  the 
leakage  field  through  the  coupling  aperture  can  be  seen  in  the  figure. 

A  second  constraint  on  the  output  cavity  is  that  it  be  stable  to  spurious  oscillations.  The 
beam  current  required  to  achieve  self-oscillation  in  the  output  cavity  is  shown  in  Fig.  2  as  a 
function  of  (uniform)  magnetic  field.  It  is  assumed  that  no  drive  power  is  applied  to  the  input 
cavity.  Curves  are  plotted  for  the  TEqh  mode  and  the  TEm  mode.  The  resonant  frequency  of  the 
fundamental  TE,,,  mode  is  near  the  cyclotron  frequency  at  4.87  GHz.  The  quality  factor  of  the 
TEjii  mode  is  72.  All  other  modes  evaluated  are  completely  stable  to  currents  below  10  A  in  the 
magnetic  field  range  indicated.  The  cavity  is  expected  to  be  unstable  to  the  operating  mode  but 
stable  to  the  fundamental  mode  at  the  design  field  of  1.83  kG.  However,  the  optimal  taper  yields 
a  field  at  the  exit  of  the  output  cavity  of  1.775  kG.  At  this  field  magnitude  the  fundamental  mode 
is  unstable  and  the  operating  mode  is  stable.  By  using  the  efficiency  code,  it  can  be  shown  that 
with  the  tapered  field  profile  the  desired  mode  is  again  unstable  and  the  TEm  mode  is  stable. 

The  evolution  of  electron  bunching  in  the  device  is  illustrated  in  Fig.  3,  where  the  distribu¬ 
tion  of  representative  particles  relative  to  the  phase  of  the  output  cavity’s  electric  field  is  plotted 
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at  various  axial  locations.  Figure  3a  represents  the  distribution  just  prior  to  the  output  cavity 
when  no  drive  power  is  applied  and  illustrates  the  initially  uniform  distribution  of  particles  in 
phase  space.  The  three  distinct  bands  of  particles  represent  the  effect  of  finite  beam  thiekness  and 
the  velocity  spread  that  results  from  canonical  angular  momentum.  Figure  3b  shows  the  distribu¬ 
tion  at  the  output  eavity  entrance  when  a  signal  of  800  W  is  injected  into  the  input  cavity.  The 
increase  in  perpendicular  velocity  spread  can  be  seen  clearly  and  is  comparable  (but  less  than) 
the  canonical  angular  momentum  spread.  The  phase  distribution  at  the  exit  of  the  output  cavity 
which  corresponds  to  the  800  W  input  case  is  shown  in  Fig  3c.  The  phase  coherence  of  the  beam 
is  essentially  destroyed  and  a  large  number  of  particles  have  given  up  over  80%  of  their  perpen¬ 
dicular  energy.  A  much  smaller  number  of  particles  have  remained  at  about  the  same  energy 
level  or  had  their  energy  inerease  slightly.  The  beam  evolution  throughout  the  cavity  is  indicated 
in  Fig.  la  where  the  r-z  projection  of  representative  particles  are  plotted.  The  figure  reveals  that 
while  some  particles  gain  energy  and  increase  their  Larmor  radii,  the  majority  of  the  particles 
lose  a  significant  amount  of  energy. 

The  simulated  results  for  the  parameters  in  Table  I  near  the  optimal  drive  power  are  summa¬ 
rized  in  Table  II.  The  nominal  electric  field  at  the  beam  is  about  8  kV/cm  in  the  input  cavity  and 
about  5  times  that  in  the  output  cavity.  The  peak  electric  field  at  the  wall  is  17  kV/cm  in  the  input 
cavity  and  zero  in  the  output  cavity.  The  total  axial  velocity  spread  includes  a  contribution  of 
about  8.5%  from  energy  spread,  which  is  reasonably  close  to  the  estimate  given  in  Eq.  (4).  The 
contribution  to  the  perpendicular  velocity  spread  from  the  energy  spread  is  about  0.8%.  The 
nominal  saturated  gain  is  about  20  dB.  The  peak  efficiency  of  53%  greatly  exceeds  the  perform¬ 
ance  of  any  second  harmonic  small-orbit  amplifier  to  date. 

The  dependence  of  efficiency  on  output  cavity  Q  is  revealed  in  Fig.  4  for  a  thin  beam  and  a 
drive  power  of  800  W.  The  output  efficiency  rises  dramatically  from  ~5%  to  -48%  as  the  quality 
factor  is  increased  from  300  to  1300.  After  that  point,  the  efficiency  only  increases  slowly  with  Q 
until  it  reaches  its  maximum  value  of  54.4  %  when  Q=2550.  The  nominal  Q  of  1825  was  se¬ 


ll 


lected  rather  than  the  optimal  value  because  its  improved  stability  and  wall  loss  properties  come 
at  only  a  small  sacrifice  in  efficiency. 

Two  drive  curves  for  parameters  near  the  optimal  configuration  (and  a  zero  thickness  beam) 
are  shown  in  Fig.  5.  The  solid  line  shows  the  peak  power  at  the  optimal  electric  field  strength. 
This  represents  the  maximum  possible  performance  but  also  requires  that  the  quality  factor  be 
varied  at  each  point.  The  zero  drive  efficiency  is  over  22%  and  requires  an  output  quality  factor 
of  3180.  The  second  curve  gives  the  output  power  for  the  original  output  cavity  with  a  Q  of  1825. 
Performance  is  poor  at  low  drive  powers,  where  the  zero  input  efficiency  is  less  than  2%.  How¬ 
ever,  the  two  curves  are  quite  similar  near  the  maximum  output  power,  which  occurs  for  the 
latter  curve  at  an  input  power  of  1.4  kW.  The  corresponding  power  and  efficiency  are  78  kW  and 
55.7%  efficiency,  respectively. 

The  dependence  of  output  efficiency  on  the  thin  beam’s  average  velocity  ratio  is  given  in 
Fig.  6.  This  ratio  is  modified  solely  by  varying  the  magnetic  field  in  the  drift  region.  The  phase 
and  amplitude  output  cavity’s  electric  field  are  optimized  at  each  point.  This  again  implies  that 
the  quality  factor  is  adjusted  to  maximize  output  power.  All  other  parameters,  however,  are  held 
at  their  nominal  values.  The  efficiency  drops  off  fairly  rapidly  with  decreasing  velocity  ratio, 
going  to  a  level  of  about  23%  at  a  velocity  ratio  near  one.  This  result  requires  a  quality  factor  of 
8570.  Efficiency  only  slowly  increases  in  the  opposite  direction,  approaching  a  level  of  57%  with 
a  velocity  ratio  of  three  and  a  quality  factor  of  1635.  This  result  would  undoubtedly  degrade  with 
the  introduction  of  finite  beam  thickness.  On  the  other  hand,  performance  could  possibly  be  im¬ 
proved  further  if  additional  optimizations  (e.g.  cavity  length  adjustment)  were  performed. 

The  effect  of  beam  thickness  on  efficiency  at  the  nominal  operating  point  has  been  examined 
in  Fig.  7.  The  efficiency  is  fairly  insensitive  to  beam  thickness  and  the  resulting  velocity  spread. 
For  example,  efficiency  barely  dips  below  50%  at  a  beam  thickness  of  0.26  cm  and  an  axial  ve¬ 
locity  spread  near  25%.  Particles  begin  to  hit  the  output  cavity  wall  when  the  thickness  exceeds 
0.28  cm. 
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The  results  from  the  single  particle  code  are  confirmed  using  the  gyroklystron  code.  In  these 
simulations,  the  pre-bunched  beam  shown  in  Fig.  3(b)  and  field  profiles  from  the  scattering  ma¬ 
trix  code  are  used  as  the  input.  The  effect  of  magnetic  field  taper  is  presented  in  Fig.  8a.  The  field 
amplitude  and  phase  are  adjusted  at  each  value  of  magnetic  taper  to  optimize  the  efficiency.  The 
beam  current  is  adjusted  such  that  the  required  quality  factor  from  the  large  signal  code  matches 
with  the  Q  value  from  the  scattering  matrix  code.  It  is  important  to  note  that  the  product  of  the 
beam  current  and  Q  is  proportional  to  the  square  of  field  amplitude  and  inversely  proportional  to 
the  efficiency.  If  the  beam  current  is  not  within  the  acceptable  range  then  the  amplitude  and 
phase  should  be  readjusted  or  the  cavity  should  be  redesigned  with  a  different  value  of  Q.  The 
efficiency  values  plotted  in  Fig.  8a,  shown  by  the  circles,  correspond  to  10  %  rms-spread  in  the 
axial  velocity.  The  corresponding  values  of  optimum  beam  current  are  shown  by  the  triangles.  It 
can  be  seen  that  the  efficiency  is  maximum  for  the  magnetic  taper  value  of  -5  G/cm  which  is 
equivalent  to  a  -3  %  taper.  The  beam  current  of  3.8  A  matches  well  with  design  current  of  4.0  A 
based  on  the  single  particle  code.  The  magnetic  taper  improves  the  efficiency  from  44.4  %  at 
zero  taper  to  5 1 .4  %.  Higher  magnetic  tapers  correspond  to  lower  values  of  field  amplitudes  re¬ 
quired  to  achieve  maximum  efficiency  and  hence  lower  values  of  optimum  beam  current.  The 
optimum  beam  current  (for  a  fixed  value  of  Q=1825)  varies  from  2.3  Amp  to  6.1  Amp  for  the 
magnetic  taper  of  -8  G/cm  to  0  G/cm.  The  corresponding  output  power  level  varies  from  40  kW 
to  95  kW. 

The  effect  of  axial  velocity  spread  on  the  design  point  is  plotted  in  Fig.  8b.  This  assumes 
that  the  beam  is  monoenergetic  at  the  entrance  of  the  output  cavity.  It  can  be  seen  from  the  fig¬ 
ure  that  the  output  efficiency  is  above  50  %  for  axial  velocity  rms-spreads  up  to  15  %.  The 
decrease  in  efficiency  with  the  velocity  spread  is  very  slow  for  this  device  as  compared  to  second 
harmonic  gyrotrons  or  other  devices.  It  is  important  to  note  that  the  velocity  spread  in  the  gyrok¬ 
lystron  code  is  “placed”  on  the  nominal  beam  distribution  and  thus  is  quite  different  from  the 
spread  computed  by  the  two  single  particle  codes.  Nonetheless,  the  agreement  between  the  two 
types  of  simulations  is  excellent. 


13 


The  final  simulation  that  we  present  utilizes  an  input  cavity  that  is  resonant  in  the  TMqio 
mode  at  half  the  output  frequency.  This  trivially  satisfies  the  requirement  in  Eq.  (1).  The  drift 
region  length  is  doubled  to  conform  to  Eq.  (2).  The  drive  power  is  adjusted  to  introduce  a  similar 
energy  spread  into  the  beam.  Most  of  the  other  parameters  remain  the  same.  All  of  the  values 
that  are  adjusted  are  summarized  in  Table  III;  those  not  adjusted  remain  as  in  Table  I.  For  a 
quality  factor  of  about  2540,  an  electronic  efficiency  of  50.9  %  is  achieved.  This  represents  only 
a  slight  decrease  (~2%)  over  the  second  harmonic  drive  case.  Furthermore,  the  reduced  input 
power  requirement,  coupled  with  the  higher  quality  factor  and  the  lower  stored  energy  boosts  the 
large-signal  gain  to  over  28  dB. 
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IV.  DISCUSSION 


The  novel  concept  of  an  axially  pre-bunched,  low-voltage,  second-harmonic,  small-orbit  gy- 
rotron  presented  here  promises  reasonable  gains  and  unprecedented  efficiencies  in  a  compact 
design.  Because  the  beam  is  linearly  streaming  throughout  most  of  the  tube,  this  device  is  easily 
more  stable  than  a  comparable  gyroklystron.  The  ability  to  use  overmoded  output  cavities  should 
enable  these  tubes  to  achieve  higher  powers  than  comparable  klystrons. 

The  specific  designs  investigated  here  utilized  a  35  kV,  4  A  beam  to  produce  over  70  kW  of 
power  at  9.9  GHz  via  an  axially  coupled  TEqi  i  output  cavity.  Simulations  were  carried  out  with  a 
single  particle  code  which  was  modified  to  accept  realistic  cold-cavity  field  profiles  from  a  scat¬ 
tering  matrix  code.  The  two-cavity  gains  with  second  (TM020)  and  first  harmonic  (TMqio)  input 
cavities  were  20  dB  and  28  dB,  respectively.  The  simulated  efficiency  for  the  X-Band  input  cav¬ 
ity  system  of  53%  included  non-ideal  effects  due  to  finite  beam  thickness  and  cusp-width.  The 
efficiency  of  the  first  harmonic  input  system  was  only  about  two  percent  lower  and  required  an 
output  cavity  Q  about  39  %  higher  than  the  previous  design.  The  addition  of  one  or  more  second- 
harmonic  buncher  cavities  should  allow  for  shorter  drift  regions,  higher  gain,  and  possibly  higher 
efficiency.  It  should  also  minimize  the  performance  discrepancies  between  the  two  input 
schemes. 

One  of  the  key  differences  between  the  axial  bunching  schemes  for  large-  and  small-orbit 
devices  is  the  axial  velocity  spread  resulting  from  the  energy  modulation.  While  the  approximate 
expression  for  velocity  spread  (Eq.  (4))  is  the  same  for  both  cases,  the  effect  of  the  input  cavity’s 
magnetic  field  appears  to  reduce  the  axial  spread  in  large-orbit  systems^'*  but  does  little  for  the 
spread  in  small  orbit  systems.  This  difference  stems  in  part  from  the  effect  of  each  electron’s  ra¬ 
dial  displacement  on  the  post  cusp  velocity  ratio.  Fortunately,  this  spread  can  often  be  minimized 
by  choosing  a  drift  length  that  is  proportional  to  the  axial  periodicity  of  the  beam’s  radial  pertur¬ 
bation.  It  is  also  possible  that  the  beam  focusing  of  the  electron  gun  could  be  designed  to  reduce 
the  velocity  spread.  Additional  work  must  be  done  on  this  concept  to  determine  its  viability. 
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Another  concept  which  requires  further  investigation  is  the  operation  of  pre-bunched  small- 
orbit  devices  above  the  second  harmonic.  Third  harmonic  operation  was  briefly  simulated  during 
this  investigation  with  a  (third-harmonic)  TM020  input  cavity  and  a  TE021  output  cavity,  but  no 
systematic  study  was  undertaken.  By  varying  the  guiding  center  radius,  we  were  able  to  fairly 
quickly  find  operating  points  with  efficiencies  exceeding  40%,  however  the  required  output  cav¬ 
ity  quality  factors  were  comparable  to  the  resistive  Qs  of  copper  cavities.  It  is  possible  that 
operation  was  hindered  by  the  induced  axial  velocity  spread,  because  harmonic  operation  is  in¬ 
creasingly  more  sensitive  to  such  spread.  However,  it  is  also  possible  that  further  optimization 
would  have  resulted  in  high  efficiency  at  a  realistic  operating  point,  particularly  since  the  mag¬ 
netic  field  resonance  is  quite  narrow  at  low  voltage.  In  systems  that  can  tolerate  higher  beam 
voltages,  the  correspondingly  broader  magnetic  field  resonances  should  facilitate  high  harmonic 
operation. 

Because  of  the  high  quality  factors  of  both  cavities  in  this  design,  the  instantaneous  band¬ 
width  of  this  device  should  be  relatively  small.  However,  we  propose  that  this  concept  can 
readily  be  extended  to  a  broad-band  version.  Because  the  beam  is  bunched  with  a  TMon  mode, 
the  input  cavity  and  drift  tube  region  can  be  replaced  by  a  helix  (or  some  other  slow  wave) 
structure.  By  the  same  token,  the  output  cavity  can  be  replaced  by  a  smooth-walled,  traveling 
wave  section.  This  device  should  realize  bandwidths  comparable  to  traveling  wave  tubes 
(TWTs),  achieve  higher  powers  and  efficiencies  than  conventional  TWTs,"*'  and  be  more  stable 
than  gyro-TWTs.'^^  The  two  regions  will  be  naturally  isolated  from  each  other  by  the  center  plug 
of  the  iron  pole  piece  (that  is  used  to  generate  the  non-adiabatic  transition).  Isolation  should  be 
further  enhanced  because  the  field  patterns  in  each  region  are  quite  distinct. 

In  the  future,  we  will  attempt  to  investigate  the  theoretical  questions  left  unanswered  in  the 
preceding  few  paragraphs.  We  also  expect  to  work  on  an  experimental  investigation  of  this  con¬ 
cept.  We  are  in  the  process  of  building  an  experimental  test  bed  which  is  energized  by  a  50  kV, 
20  A,  5  ps  pulse-line  modulator.  We  have  the  capability  to  produce  magnetic  fields  of  about  2 
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kG  over  the  required  distance  and  are  currently  designing  an  electron  gun  compatible  with  the 
beam  parameters  of  this  design. 

This  work  was  supported  by  the  Tri-Services  Program  for  Vacuum  Electronics. 
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Table  I.  Second  Harmonic  Amplifier  System  Parameters 


Beam  Parameters 


Voltage  (kV) 

35 

Current  (A) 

4 

(before  cusp) 

(after  cusp) 

Average  radius  (cm) 

1.50 

1.18 

Beam  thickness  (mm) 

1.00 

7.33 

Velocity  ratio  {v\Jv^ 

0.00 

2.28 

Input  Cavity 

Drive  frequency  (GHz) 

9.9 

Operating  mode 

TMo20 

Resistive  Quality  factor  (Q) 

2800 

Radius  (cm) 

2.66 

Length  (cm) 

0.20 

Magnetic  Field  Parameters 

Buncher  magnetic  field  (kG) 

0.985 

Output  magnetic  field  (kG) 

1.83 

Cyclotron  frequency  (GHz) 

4.794 

Cusp  width  (cm) 

0.50 

Cavity  -  cusp  spacing  (cm) 

13 

Output  field  taper 

-3% 

Output  cavity 

Output  frequency  (GHz) 

9.9 

Operating  mode 

H 

o 

Diffractive  Q 

1825 

Resistive  Q 

28,580 

(main  section) 

(cavity  lip) 

Radius  (cm) 

1.864 

1.70 

Length  (cm) 

10.000 

1.09 
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Table  11.  Second  Harmonic  Amplifier  Performance 


Input  cavity  results 

Input  drive  power  (W) 

800 

at  beam  (kV/cm) 

7.99 

Drift  region  results 

Total  Av^  (%) 

12.2 

Total  Avx  (%) 

1.8 

Energy  spread  (%) 

3.2 

Output  cavity  results 

Power  (kW) 

74.2 

Efficiency  (%) 

53.0 

Gain  (dB) 

19.7 

E^  at  beam  (kV/cm) 

39 
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Table  III.  First  Harmonic  Drive  System  Parameters 


Input  Cavity 

Drive  frequency  (GHz) 

4.95 

Operating  mode 

TMqio 

Resistive  Quality  factor  (Q) 

3630 

Radius  (cm) 

2.318 

Length  (cm) 

0.40 

Input  drive  power  (W) 

110 

Ez  at  beam  (kV/cm) 

4.1 

Drift  region 


Cavity  -  cusp  spacing  (cm) 

26 

Total  AVz  (%) 

11.8 

Total  Avx  (%) 

1.7 

Energy  spread  (%) 

3.3 

Output  cavity 


Diffractive  Q 

2540 

Power  (kW) 

71.3 

Efficiency  (%) 

50.9 

Gain  (dB) 

28.1 

E^  at  beam  (kV/cm) 

44 
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Figure  Captions 


1.  The  axially-modulated,  cusp-injected,  small-orbit  gyrotron;  (a)  system  schematic,  (b)  axial 
magnetic  field  profile,  and  (c)  normalized  electric  field  profile  at  the  beam  radius  (E^  in  the 
input  cavity  and  in  the  output  cavity). 

2.  The  output  cavity  start-oscillation  curves.  The  squares  indicate  the  TEqh  mode  and  the  cir¬ 
cles  represent  the  TEm  mode. 

3.  Representative  phase  distribution  of  electrons:  (a)  output  cavity  entrance  (Pjn^O  W),  (b)  out¬ 
put  cavity  entrance  (Pin=800  W),  and  (c)  output  cavity  exit  (Pin=800  W). 

4.  The  dependence  of  output  efficiency  on  the  output  cavity  quality  factor  (Pin=800  W). 

5.  System  drive  curves.  The  solid  line  shows  the  output  efficiency  at  the  optimal  electric  field 
strength.  The  dashed  line  gives  the  output  efficiency  for  a  constant  quality  factor  (Q=l  825). 

6.  The  dependence  of  output  efficiency  on  average  velocity  ratio. 

7.  The  dependence  of  output  efficiency  on  beam  thickness. 

8.  Gyroklystron  code  results:  (a)  output  efficiency  (solid  line)  and  optimal  current  (dashed  line) 
as  a  function  of  magnetic  field  taper  and  (b)  output  efficiency  as  a  function  of  velocity 
spread. 


1.  The  axially-modulated,  cusp-injected,  small-orbit  gyrotron:  (a)  system  schematic,  (b)  axial 
magnetic  field  profile,  and  (c)  normalized  electric  field  profile  at  the  beam  radius  (E^  in  the 
input  cavity  and  in  the  output  cavity). 
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6.  The  dependence  of  output  efficiency  on  average  velocity  ratio. 
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7.  The  dependence  of  output  efficiency  on  beam  thickness. 
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8.  Gyroklystron  code  results:  (a)  output  efficiency  (solid  line)  and  optimal  current  (dashed  line) 

as  a  function  of  magnetic  field  taper  and  (b)  output  efficiency  as  a  fimction  of  velocity 
spread. 
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Design  of  an  Efficient,  Low  Voltage,  Third  Harmonic,  Large-Orbit  Gyrotron  Amplifier 

with  a  Vane-Resonator  Output  Cavity 
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Abstract 

We  present  the  design  of  a  100  kW,  J-band,  third  harmonic  large-orbit  gyrotron 
amplifier  which  utilizes  the  interaction  between  a  45  kV,  4  A  beam  and  a  vane  resonator 
output  cavity  operating  in  the  “7r”-mode.  An  efficiency  of  55%  is  predicted  with  a  large 
signal  gain  near  20  dB  by  a  single  particle  code  which  takes  into  account  non-ideal  ef¬ 
fects  associated  with  finite  beam  thickness  and  finite  magnetic  field  transition  widths. 
High  efficiency  is  achieved  by  velocity  modulation  of  an  axially-streaming  annular 
beam  via  a  short  TM310  drive  cavity.  Ballistically  created  axial  bunches  are  converted 
into  azimuthal  bunches  when  the  beam  encounters  a  non-adiabatic,  balanced  magnetic 
field  reversal  at  the  end  of  a  30  cm  drift  region.  The  design  of  this  tube  is  presented  and 
its  performance  is  completely  characterized  before  the  prospects  for  the  operation  of  this 
low  voltage  configuration  at  other  harmonics  are  explored. 


I.  INTRODUCTION 


Microwave  generation  at  high  harmonics  of  the  cyclotron  frequency  from  large-orbit  gyro- 
trons  has  undergone  considerable  theoretical  and  experimental  scrutiny  in  the  past  two  decades  at 
a  number  of  institutions  [1-16].  Because  of  the  reduced  magnetic  field  requirements,  these  de¬ 
vices  potentially  have  a  considerable  advantage  over  conventional  first  harmonic  gyrotrons  [17- 
19]  in  a  number  of  applications  that  require  compact,  lightweight  tubes.  The  chief  disadvantage 
of  high  harmonic  operation  is  the  reduction  in  efficiency  that  usually  accompanies  any  increase 
in  harmonic  number.  One  possible  way  to  compensate  for  this  weakness  is  to  use  energy  recov¬ 
ery  techniques  to  restore  the  net  efficiency  to  a  level  comparable  to  the  first  harmonic  devices 
[13,20,21].  Unfortunately,  this  approach  adds  to  the  complexity  of  the  power  supply  and  micro- 
wave  output  system.  Also,  there  is  typically  a  large  capital  cost  associated  with  the  collector 
hardware  (insulators,  electrodes,  magnets,  etc.).  Furthermore,  this  approach  does  nothing  to  re¬ 
duce  cathode  loading  requirements  or  space-charge  effects  in  the  gun  and  beam  transport  system. 

An  alternate  approach  to  enhanced  efficiency  involves  pre-bunching  the  beam  to  achieve  a 
non-uniform  phase-space  distribution  prior  to  the  entrance  of  the  output  cavity.  Several  such  ef¬ 
ficiency  enhancement  techniques  have  been  attempted  in  linear  devices  [22-24],  and  all  of  them 
may  be  applicable  to  those  large-orbit  devices  which  introduce  perpendicular  energy  into  the 
beam  via  a  non-adiabatic,  magnetic  field  reversal.  One  of  these  techniques  utilizes  a  modified 
version  of  the  klystron’s  ballistic  bunching  process  to  generate  the  required  velocity  distribution. 
A  device  based  on  this  technique,  which  we  call  the  axially-modulated,  cusp-injected,  large-orbit 
gyrotron,  has  been  explored  in  an  earlier  paper  for  a  system  with  a  100  kV,  25  A  beam  [25]. 

A  schematic  of  this  device  is  shown  in  Fig.  1.  The  axial  bunches  are  created  ballistically  af¬ 
ter  the  beam  passes  through  an  input  cavity  which  is  driven  with  a  TM^^iq  circularly-polarized 
mode.  The  axial  bunches  are  converted  into  azimuthal  bunches  when  the  beam  encounters  the 
non-adiabatic  magnetic  transition  (cusp).  If  the  drive  frequency  is  at  the  harmonic  of  the  cy¬ 
clotron  frequency,  the  net  effect  is  that  the  beam  will  form  m  spokes  that  will  rotate  at  the  cyclo¬ 
tron  frequency.  If  an  output  cavity  that  operates  in  a  TEn,np  mode  is  placed  immediately  after  the 
cusp,  efficient  energy  extraction  can  occur.  For  example,  the  previous  design  study  demonstrated 
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that  an  efficiency  of  40%  was  achievable  in  X-Band  at  the  fourth  harmonic  via  an  interaction 
with  a  right-circular  TE4ji  cavity. 

While  the  result  quoted  in  the  previous  paragraph  appears  promising,  there  are  some  practi¬ 
cal  limits  to  the  applicable  range  of  parameters  for  this  device.  Of  particular  concern  is  the  qual¬ 
ity  factor  required  to  produce  the  optimal  field  strength  in  a  low  voltage,  low  current  device 
operating  at  a  high  harmonic  number.  For  example,  the  radius  of  the  maximum  electric  field  in  a 
cavity,  relative  to  the  wall  radius,  increases  with  increasing  harmonic  number.  Thus,  the  higher 
the  harmonic  number  the  further  out  the  beam  should  be.  However,  the  beam  radius  of  the  large- 
orbit  configuration  is  equal  to  the  Larmor  radius,  which,  for  a  particular  cyclotron  frequency,  is 
proportional  to  the  perpendicular  velocity.  Thus,  the  beam  radius  in  low  voltage  systems  de¬ 
creases  fairly  rapidly  with  decreasing  beam  energy.  These  two  seemingly  incompatible  require¬ 
ments  often  combine  to  yield  an  unrealistic  requirement  on  the  output  cavity  quality  factor. 

There  are  at  least  two  ways  to  overcome  this  difficulty  and  design  efficient,  low  voltage, 
harmonic  devices.  One  way  is  to  use  a  pre-bunched,  small-orbit  beam.  Axial  modulation  is  a  vi¬ 
able  technique  for  this  approach  because  a  non-reversal,  non-adiabatic  magnetic  transition  can  be 
use  to  generate  the  small  orbit  beam.  The  guiding  center  radius  can  be  adjusted  to  place  the  beam 
at  a  location  where  the  electric  field  is  relatively  strong.  This  approach  has  been  investigated 
previously  via  the  design  of  a  second  harmonic  small-orbit  system  [26].  In  this  design,  a  35  kV,  4 
A  beam  was  pre-bunched  by  a  short  TM020  cavity  and  interacted  with  a  10  cm  long  TEqh  output 
cavity  to  theoretically  produce  over  70  kW  of  power  in  X-Band.  The  simulated  efficiency  and 
large-signal  gain  were  approximately  53%  and  20  dB,  respectively,  when  the  beam  thickness  and 
finite  magnetic  transition  width  were  taken  into  consideration.  At  approximately  1800,  the  re¬ 
quired  output  cavity  quality  factor  was  well  below  the  resistive  Q  of  28,600  and  the  cavity  was 
stable  to  all  spurious  modes.  However,  the  simulations  indicated  that  the  quality  factor  of  a  com¬ 
parable  third  harmonic  system  would  be  too  large. 

The  second  approach  is  to  use  a  vane-resonator  to  couple  various  azimuthal  modes  together. 
The  resultant  cavity  modes  can  support  the  necessary  azimuthal  field  variation  with  a  much 
lower  cutoff  frequency.  Thus,  the  wall  radius  can  be  reduced  considerably  as  compared  to  the 
smooth  wall  case  and  the  interaction  impedance  can  be  significantly  increased.  This  is  the  ap¬ 
proach  that  we  will  investigate  in  this  paper.  In  the  next  section  we  will  discuss  the  basic  ampli- 
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fier  configuration  and  describe  the  computer  codes  used  in  the  analysis.  In  Section  III,  we  will 
present  the  results  of  a  third  harmonic  example  which  utilizes  a  45  kV,  4  A  beam  to  efficiently 
generate  100  kW  of  power  at  5.5  GHz.  These  parameters  were  selected  to  facilitate  a  future  proof 
of-principle  experiment  in  our  laboratory.  In  addition  to  characterizing  tube  stability  and  large- 
signal  characteristics,  we  will  also  discuss  the  sensitivity  of  the  device  to  parameter  variations. 
The  results  of  this  study  and  prospects  for  high  harmonic  operation  are  summarized  in  the  final 
section. 
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II.  DEVICE  OPERATION  AND  MODELING 


A  detailed  description  of  the  operation  of  this  device  was  given  in  [25].  In  this  section  we 
briefly  note  the  key  aspects  of  the  design  and  simulation  processes  and  detail  only  the  novel  as¬ 
pects  of  this  work,  which  involve  the  use  of  a  vane  resonator  (VR)  output  cavity. 

In  order  to  maximize  the  circuit  gain,  the  beam  must  be  placed  near  the  maximum  of  the 
axial  electric  field  in  the  input  cavity,  which  is  increasingly  close  to  the  wall  radius  with  increas¬ 
ing  harmonic  number.  The  wall  radius  in  a  TM^jq  cavity  is  proportional  to  the  drive  signal’s 
wavelength  while  the  beam  radius  is  tied  to  the  magnetic  field  and  the  beam  energy.  This  prob¬ 
lem  is  easier  to  remedy  than  the  analogous  one  with  the  output  cavity:  either  the  input  cavity  can 
be  moved  into  the  magnetic  field’s  adiabatic  compression  region  where  the  beam  radius  is  larger 
(if  one  exists)  or  dielectrics  can  be  inserted  into  the  input  cavity  to  decrease  the  wall  radius.  We 
have  selected  the  latter  approach  for  this  design.  Because  the  cavity  is  short  and  the  beam  is  thin 
and  streaming  linearly  in  that  region,  there  should  be  no  problem  with  charge  build-up.  We 
model  the  cavity  with  a  code  that  assumes  three  dielectric  regions  radially  and  calculates  the 
quality  factor  assuming  a  good  conducting  (copper)  wall  and  an  ideal  dielectric. 

The  required  drift  tube  length  d  depends  on  the  magnitude  of  the  energy  spread  (AE/Eg)  im¬ 
parted  to  the  beam: 

d=^^X/i2AE/E,)  (1) 

where  Po  is  the  streaming  velocity  (normalized  to  the  speed  of  light)  and  X  is  the  EM  wavelength 
(of  the  input  cavity  signal).  The  magnetic  transition  at  the  end  of  the  drift  region  is  shortened  by 
a  thin  iron  pole  piece.  Long  solenoids  on  either  side  of  the  pole  piece  generate  approximately 
constant  axial  magnetic  fields  in  each  region  (Bg  on  the  input  side  and  Bo  on  the  output  side).  The 
axial  magnetic  field  is  modeled  by  a  piece-wise  linear  function  and  the  radial  field  is  given  by 
Maxwell’s  equations.  The  spread  in  beam  energy  translates  into  additional  velocity  spread  after 
the  magnetic  transition.  When  the  effects  of  the  input  cavity’s  magnetic  fields  are  neglected,  the 
induced  spread  can  be  approximated  by  [25]: 
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In  (2),  Yo  is  the  relativistic  mass  factor. 

A  single  particle  code  is  used  to  simulate  the  beam’s  motion  through  the  input  cavity,  the 
drift  region,  and  the  non-adiabatic  magnetic  transition.  Given  the  beam  parameters,  electrons  are 
initially  “launched”  over  a  range  of  axial  positions,  times,  and  angles  that  represent  the  entirety 
of  phase  space  (after  taking  into  account  any  symmetries  in  the  system).  The  initial  velocities  of 
the  particles  are  assumed  to  be  purely  axial.  A  minimum  of  100  particles  are  used  to  model  the 
beam  when  zero  radial  beam  thickness  is  assumed  and  more  than  300  particles  are  used  to  model 
finite  thickness  beams.  After  a  Runge-Kutta  start-up,  a  predictor-corrector  integration  scheme  is 
used  to  solve  the  equations  of  motion.  Convergence  is  checked  by  varying  the  time  step  and  the 
number  of  particles  and  by  verifying  conservation  of  energy  and  canonical  angular  momentum 
for  each  particle.  This  “buncher”  code  produces  the  phase  space  distributions  of  the  electrons  at 
the  end  of  the  transition  region,  calculates  the  key  beam  statistics  {e.g.  the  average  velocity  ratio 
and  guiding  center  radius  and  the  spreads  in  Vj^,  v^.,  and  energy),  and  indicates  the  effectiveness 
of  the  bunching  process. 

A  cross-section  of  the  output  cavity  is  displayed  in  Fig.  2.  The  wall  radius  is  denoted  r^  and 
the  slot  depth  is  dj  so  that  the  slot  radius  is  rj  =  r^  +  dj.  We  assume  that  there  are  slots,  so  that 
the  angular  span  of  one  period  is  Qf  =  2nl  n^.  The  angular  opening  of  each  slot  is  A0.  To  solve  for 
the  dispersion  relation  and  subsequently  the  cutoff  wavenumbers  for  the  VR  structure,  we  make 
the  usual  approximation  of  a  constant  electric  field  across  the  slot  opening  and  match  the  mag- 
netie  field  there  only  on  average  [13].  If  we  let  the  phase  advance  between  adjacent  slots  be: 

(p^=2nljn^  (3) 

where  is  an  integer  in  the  range  ,  then  the  dispersion  relation  is  given  by: 
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Here  and  Xg  -  for  ^  =  ^(o)  /  c)  -k]  ^  and  where  co  is  the  angular  frequency,  c  is  the 

speed  of  light,  and  is  the  axial  wavenumber.  Jn(x)  and  Yn(x)  are  Bessel  functions  of  the  first 
and  second  kind,  respectively. 

A  simple  code  is  used  to  find  the  cutoff  frequency  of  a  VR  structure  from  (4),  given  the  di¬ 
mensions  listed  above  and  a  particular  value  of  4.  The  usual  two  operating  modes  are  the  “n”- 
mode  {Ig  =  nj  2)  and  the  “27T”-mode  (/„  =  0).  We  will  only  consider  the  “7:”-mode  in  this  paper 
because  it  has  in  general  a  lower  start-oscillation  threshold  than  the  other  mode.  The  EM  fields 
are  generated  by  specifying  an  axial  length  for  the  cavity  and  assuming  that  the  axial  variation  is 
sinusoidal  (the  closed  cavity  approximation).  The  resistive  quality  factor  is  calculated  by  assum¬ 
ing  a  good  conductor  (copper)  wall,  evaluating  the  magnetic  field  at  all  surfaces,  and  integrating 
the  wall  losses  either  analytically  or  numerically.  The  actual  resistive  Q  will  presumably  be 
somewhat  larger,  since  the  endwalls  will  be  partially  removed. 

A  second  single  particle  code  takes  the  results  from  the  “buncher”  code  and  uses  the  same 
numerical  techniques  to  calculate  the  extraction  efficiency  of  the  EM  wave  in  the  output  cavity 
under  steady-state  conditions.  The  “efficiency”  code  also  takes  the  results  from  the  output  cavity 
code  to  generate  the  electromagnetic  fields.  The  applied  dc  magnetic  field  in  the  output  cavity  is 
assumed  to  vary  linearly.  This  “efficiency”  code  is  iterated,  while  varying  the  cavity  and  mag¬ 
netic  field  parameters,  imtil  an  optimal  design  is  achieved.  During  each  run,  the  amplitude  and 
phase  of  the  wave  are  adjusted  until  maximum  efficiency  is  obtained.  In  addition  to  the  effi¬ 
ciency,  the  code  reports  the  required  diffractive  Q.  By  decreasing  the  field  amplitude  in  the  out¬ 
put  cavity  towards  zero,  this  code  is  used  to  calculate  the  start  currents  for  the  various  modes  in 
tapered  magnetic  fields.  The  code  has  been  benchmarked  successfully  against  the  start-oscillation 
results  in  [8]  by  using  an  unbunched  beam. 
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III.  SIMULATION  RESULTS  FOR  A  THIRD  HARMONIC  DESIGN 


The  parameters  for  a  third  harmonic  large-orbit  design  are  given  in  Table  I.  The  system  is 
designed  to  produce  microwaves  at  5.5  GHz  via  the  interaction  between  a  45  kV,  4  A  beam  and 
the  microwave  circuit.  Given  the  beam  energy  and  the  expected  magnetic  field  level,  the  beam 
radius  of  1.125  cm  is  required  to  provide  a  velocity  ratio  near  2.2.  This  beam  can  be  generated 
with  a  thin  annular  cathode  and  a  modest  magnetic  compression.  For  example,  a  compression  of 
about  seven  would  result  in  a  cathode  loading  of  slightly  above  1  AJcm,  an  average  emitter  ra¬ 
dius  of  3  cm,  and  a  cathode  thickness  of  2  mm.  The  radial  thickness  indicated  corresponds  to  an 
ideal  post-cusp  velocity  spread  from  canonical  angular  momentum  considerations  of  about  Avj^= 
1.9%.  The  dc  space  charge  depression  in  the  entire  circuit  is  negligible. 

The  drive  cavity  is  a  dielectric-loaded,  right-circular  TM310  cavity  with  the  dimensions  indi¬ 
cated  in  Table  I.  The  cavity  width  is  selected  to  give  a  beam  coupling  factor  [27]  of  about  0.98.  A 
dielectric  with  a  relative  permittivity  of  about  13  is  assumed  to  be  placed  above  and  below  the 
beam  location,  leaving  a  clearance  of  at  least  1  mm.  The  cavity  radius  is  adjusted  to  achieve  the 
desired  operating  frequency.  With  these  parameters,  the  axial  electric  field  shown  in  Fig.  3  is 
realized.  Thus,  the  beam  is  located  in  a  region  where  the  electric  field  is  nearly  constant  and  op¬ 
timal.  The  azimuthal  magnetic  field  is  also  indicated  in  the  figure.  This  field  is  quite  small  near 
the  beam  radius  and  has  little  effect  on  particle  motion.  The  radial  magnetic  field  is  similar  in 
profile  to  the  electric  field  and  is  not  shown.  There  are  a  number  of  good  dielectrics  with  relative 
permittivities  near  10  £ind  we  expect  that  a  real  device  would  use  one  of  these  and  would  also 
possibly  place  the  cavity  near  the  end  of  the  magnetic  compression  region  to  maximize  the  gain. 

The  dc  magnetic  field  profile  throughout  the  circuit  is  determined  during  the  efficiency  op¬ 
timization  procedure.  The  axial  magnetic  field  at  the  entrance  to  the  output  cavity  is  such  that 
three  times  the  cyclotron  frequency  is  about  14%  lower  than  the  drive  frequency.  The  13%  upta- 
per  leaves  the  drive  frequency  nearly  equal  to  thrice  the  cyclotron  frequency  at  the  output  cavity 
end.  The  magnetic  transition  is  not  perfectly  balanced;  the  field  strength  in  the  drift  region  is 
about  8%  lower  than  the  initial  field  after  the  transition.  The  distance  between  the  input  cavity 
and  cusp  transition  is  optimized  with  the  numerical  codes,  but  the  initial  length  is  chosen  from 
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(1)  so  that  maximal  bunching  could  be  achieved  with  an  energy  spread  under  4%.  The  transition 
width  should  readily  achievable  with  the  aid  of  one  or  more  iron  pole  pieces. 

The  length  of  the  output  cavity  is  varied  to  maximize  efficiency  while  keeping  the  required 
Q  less  than  1 0%  of  the  resistive  Q  of  a  copper  cavity.  This  restriction  allows  for  extraction  of 
about  90%  of  the  microwave  power.  The  chosen  length  is  such  that  unperturbed  particles  un¬ 
dergo  approximately  414  revolutions  inside  the  cavity.  The  cavity  radius  is  selected  to  minimize 
beam  interception  and  the  slot  depth  is  adjusted  to  produce  the  required  resonant  frequency.  Six 
vanes  are  required  for  “7r”-mode  operation.  The  slot  opening  is  fixed  at  50%  of  the  angular  pe¬ 
riod.  The  azimuthal  electric  field  and  axial  magnetic  field  radial  profiles  are  plotted  in  Fig.  4  at  a 
fixed  azimuthal  location  which  corresponds  to  the  center  of  a  vane.  They  are,  of  course,  out  of 
time  phase  by  90°.  The  discontinuity  in  the  curves  at  r^^  results  from  the  truncation  of  the  infinite 
series  in  (4).  The  azimuthal  electric  field  at  the  initial  average  beam  radius  is  nearly  half  of  the 
peak  field. 

The  beam  current  required  to  achieve  self-oscillation  in  the  “Tc”-mode  of  the  output  cavity  is 
shown  in  Fig.  5  as  a  function  of  (uniform)  magnetic  field.  It  is  assumed  that  no  drive  power  is 
applied  to  the  input  cavity.  The  resonant  frequencies  of  the  various  azimuthal  modes  are  given  in 
Table  II.  The  “'n:”-mode  is  the  only  unstable  mode  in  the  magnetic  field  range  from  600  to  700  G. 
However,  even  the  “7t”-mode  is  stable  for  the  optimal  tapered  magnetic  field  profile  at  the  de¬ 
sign  current  of  4  A.  The  optimal  zero-drive  efficiency  is  found  to  be  27%  at  a  starting  magnetic 
field  of  662  G  and  a  weak  uptaper.. 

The  simulated  results  for  the  parameters  in  Table  I  near  the  optimal  drive  power  are  summa¬ 
rized  in  Table  III.  The  nominal  electric  field  at  the  average  beam  location  is  about  1 1  kV/cm  in 
the  input  cavity  and  almost  double  that  in  the  output  cavity.  The  total  axial  velocity  spread  in¬ 
cludes  a  contribution  of  about  1.9%  from  energy  spread,  which  is  substantially  below  the  esti¬ 
mate  given  in  (2)  because  of  ac  magnetic  field  effects.  The  contribution  to  the  perpendicular 
velocity  spread  from  the  energy  spread  is  about  1 .7%.  The  bunching  effectiveness  is  quantified 
by  computing  the  RMS  angular  distribution  of  particles  at  the  entrance  to  the  output  cavity,  rela¬ 
tive  to  the  phase  of  the  EM  wave.  This  value  is  scaled  and  subtracted  from  1 00%  to  yield  a  value 
of  0%  for  a  uniform  beam  and  a  value  of  100%  for  a  perfectly  bunched  beam.  The  nominal  satu¬ 
rated  gain  is  just  above  19  dB.  The  simulated  efficiency  of  55%  is  more  than  twice  the  maximum 
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unbunched  efficiency  and  is  comparable  to  or  better  than  the  optimal  theoretical  efficiencies  of 
most  first  harmonic  gyrotrons.  With  a  180  kW  beam,  about  100  kW  of  microwave  power  should 
be  produced. 

The  axial  evolution  of  electron  bunching  in  the  device  at  the  optimal  point  is  illustrated  in 
Fig.  6,  where  the  r-z  projection  of  the  trajectories  of  representative  particles  are  plotted.  The  in¬ 
put  cavity  entrance  is  at  z=0.  The  scalloping  of  the  beam  in  the  drift  region  is  due  to  the  Lorentz 
force  brought  about  by  the  radial  magnetic  field  in  the  input  cavity.  Since  the  guiding  center  ra¬ 
dius  of  the  beam  is  nearly  zero,  a  particle’s  radial  location  is  a  good  indication  of  its  energy.  The 
interaction  between  the  wave  and  the  beam  is  quite  strong,  with  many  particles  gaining  energy  in 
the  first  half  of  the  cavity  only  to  lose  it  in  the  second  half.  Clearly  the  cavity  radius  is  just  large 
enough  to  avoid  beam  scrape-off  A  larger  radius  could  be  used,  but  the  required  quality  factor 
would  increase. 

Figure  7a  shows  the  distribution  of  representative  particles  in  the  transverse  plane  (for  the 
same  optimal  parameters  of  Fig.  6)  at  the  entrance  to  the  output  cavity.  Though  it  is  difficult  to 
quantify  the  angular  density  distribution,  three  distinct  particle  bunches  are  clearly  visible.  The 
resultant  phase  distribution  at  the  exit  of  the  output  cavity  is  shown  in  Fig  7b.  The  phase  coher¬ 
ence  of  the  beam  is  essentially  destroyed  and  a  large  number  of  particles  have  given  up  the  ma¬ 
jority  of  their  perpendicular  energy.  A  much  smaller  number  of  particles  have  remained  at  about 
the  same  energy  level  or  had  their  energy  increase  slightly. 

The  efficiency  is  plotted  in  Fig.  8  as  a  function  of  beam  thickness.  Also  indicated  in  the  fig¬ 
ure  is  the  axial  velocity  spread  computed  by  the  “buncher”  code.  The  maximum  efficiency  of 
over  67%  is  obtained  with  a  zero  thickness  beam.  The  efficiency  drops  steadily  to  about  55%  at 
the  nominal  design  thickness  and  begins  to  fall  off  more  rapidly  after  that  point.  The  bunching 
effectiveness  (not  plotted)  is  essentially  independent  of  beam  thickness.  The  required  quality 
factor  varies  slowly  from  1070  at  zero  thickness  to  1240  at  the  0.75  mm  point  and  then  jumps  to 
1 520  at  the  final  point.  The  increase  in  velocity  spread  (from  canonical  angular  momentum  con¬ 
servation)  is  believed  to  be  the  main  cause  for  the  decrease  in  efficiency.  Beam  scrape-off  is 
probably  also  a  contributing  factor  to  the  efficiency  degradation,  particularly  at  the  final  thick¬ 
ness  point.  If  a  shielded  electron  source,  such  as  the  Advanced  Center-Post  gun  [28],  were  used 
instead  of  an  immersed  emitter,  efficiencies  for  a  finite  thickness  beam  could  approach  67%.  In 
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the  remainder  of  this  section  we  present  mainly  the  characteristics  of  the  zero-thickness  (or 
shielded)  beam;  the  characteristics  of  the  0.75  mm  beam  are  found  to  be  similar,  with  a  reduction 
in  efficiency  of  about  10%  near  the  optimal  power  points. 

The  drive  curve  near  the  optimal  system  parameters  (but  with  a  zero  thickness  beam)  is 
shown  in  Fig.  9.  The  relative  bunching  effectiveness  is  given  in  the  figure  by  the  dashed  line.  The 
curve  is  generated  assuming  a  constant  quality  factor  in  the  output  cavity  of  about  1100.  The  de¬ 
vice  is  zero-drive  stable  and  enjoys  a  fairly  broad  saturation  region.  The  induced  energy  spread  in 
the  beam  is  about  6.2%  at  the  maximum  drive  power  of  3  kW. 

The  dependence  of  efficiency  on  output  cavity  Q  is  displayed  in  Fig.  10.  All  other  parame¬ 
ters  are  held  fixed  at  the  optimal  system  configuration  values.  The  output  efficiency  rises  steadily 
from  zero  to  about  45%  as  the  quality  factor  is  increased  to  875.  After  that  point,  the  efficiency 
increases  more  rapidly  with  Q  until  it  reaches  its  maximum  value  of -1100.  Electric  fields  that 
are  indicative  of  quality  factors  above  1150  result  in  beam  scrape-off  and  consequently  are  not 
plotted.  Undoubtedly,  scrape-off  could  be  avoided  and  comparable  efficiencies  achieved  if  we 
were  to  adjust  some  of  the  other  parameters  (e.g.  magnetic  field  taper). 

The  dependence  of  output  efficiency  on  the  thin  beam’s  average  velocity  ratio  is  given  in 
Fig.  11.  This  ratio  is  modified  solely  by  varying  the  magnetic  field  in  the  drift  region;  the  re¬ 
quired  field  at  each  point  is  indicated  in  the  figure  by  the  dashed  line.  The  phase  and  amplitude 
of  the  output  cavity’s  electric  field  are  optimized  at  each  point.  This  implies  that  the  quality  fac¬ 
tor  is  adjusted  to  maximize  output  power.  All  other  parameters,  however,  are  held  at  their  nomi¬ 
nal  values.  The  first  maximum  in  efficiency  is  67%,  occurs  near  the  nominal  velocity  ratio  of  2.2, 
and  drops  off  fairly  rapidly  with  decreasing  velocity  ratio,  going  to  a  level  of  about  13%  at  a  ve¬ 
locity  ratio  near  one.  The  maximum  quality  factor  of  1850  is  required  at  the  minimum  efficiency 
point.  There  is  a  dip  in  the  efficiency  near  a  velocity  ratio  2.35  before  it  returns  to  the  67%  level 
at  a  velocity  ratio  of  2.5.  A  buncher  field  of  610  Gauss  yields  a  velocity  ratio  above  2.8  and  re¬ 
sults  in  particle  reflections  in  the  output  cavity’s  tapered  magnetic  field. 
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IV.  DISCUSSION 


In  this  paper  we  have  demonstrated,  via  a  concrete  design,  the  feasibility  of  the  low  voltage, 
high  efficiency,  pre-bunched,  gyrotron  concept  at  the  third  harmonic.  The  principle  extension  of 
this  work  beyond  our  previous  study  [25]  involves  the  introduction  of  a  vane-resonator  output 
cavity.  This  configuration  allows  a  reduction  in  the  transverse  dimensions  of  the  cavity  for  a 
given  frequency  and  azimuthal  mode  number  and  results  in  vastly  improved  beam  coupling  to 
low  energy,  axis-encircling  beams.  With  this  approach,  we  are  able  to  produce  simulated  effi¬ 
ciencies  of  67%  with  a  shielded  beam  and  55%  with  an  immersed  beam  (which  had  an  axial  ve¬ 
locity  spread  of  about  9%).  With  a  beam  voltage  of  45  kV  and  a  current  of  4  A,  the  resultant 
power  is  100  kW.  The  design  is  most  sensitive  to  the  required  output  cavity  quality  factor,  which 
at  the  optimal  point  is  acceptable  at  less  than  10%  of  the  output  cavity’s  resistive  Q. 

We  have  also  taken  a  preliminary  look  at  efficiency  enhancement  in  higher  harmonic  systems 
with  comparable  beam  voltages.  In  Fig.  12a,  we  plot  the  simulated  efficiency  of  the  large-orbit 
device  as  a  function  of  harmonic  number.  The  solid  line  indicates  the  simulated  efficiency 
achievable  for  a  zero-thickness  beam  and  the  dashed  line  gives  the  0.75  mm  beam  result.  The 
third  harmonic  point  corresponds  to  the  system  described  in  the  previous  section.  For  the  other 
harmonic  results,  the  magnetic  field  and  beam  parameters  are  fixed  at  the  values  given  in  Table  1. 
However,  the  cavity-cusp  spacing,  and  the  cavities’  harmonic  number,  radial  dimensions,  and 
number  of  (output)  vanes  are  adjusted  to  the  appropriate  values.  Furthermore,  the  output  cavity 
quality  factor  is  adjusted  to  optimize  efficiency.  With  the  exception  of  the  second  harmonic  re¬ 
sult,  the  finite  beam  thickness  results  in  an  efficiency  decrease  of  about  10  -  12%. 

In  Fig.  12b  we  plot  the  required  output  quality  factors  for  the  thin  (solid  line)  and  thick 
(dashed  line)  beams.  Operation  above  the  fourth  harmonic  is  not  practical  because  the  required  Q 
is  comparable  to  or  exceeds  the  resistive  Q  of  a  copper  cavity  (given  by  the  dot-dashed  line  in  the 
figure).  These  results  for  the  higher  harmonics  can  undoubtedly  be  improved  upon  by  a  more 
complete  optimization  (e.g.  adjustment  of  the  output  cavity  length),  but  the  basic  trend  would 
still  persist.  Thus,  operation  at  harmonics  significantly  above  the  fourth  will  require  an  increase 
in  beam  voltage  and  /  or  power. 

In  the  near  future,  we  will  attempt  a  proof-of-principle  experiment  based  on  this  concept.  We 
are  in  the  process  of  building  a  test  bed  which  is  energized  by  a  50  kV,  20  A,  5  ps  pulse-line 
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modulator.  The  frequency  in  this  study  is  selected  to  be  compatible  with  sources  and  magnetic 
field  coils  which  are  readily  available  to  us.  We  are  currently  designing  an  electron  gun  that  is 
compatible  with  the  beam  parameters  of  this  design. 

This  work  was  supported  by  the  Department  of  Defense  Tri-Services  Program  for  Vacuum 
Electronics. 
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Table  I.  Third  Harmonic  Amplifier  System  Parameters. 


Beam  Parameters 

Voltage  (kV) 

45 

Current  (A) 

4 

Average  radius  (cm) 

1.125 

Beam  thickness  (mm) 

0.75 

Streaming  velocity 

0.394 

Velocity  ratio  (v_l/v^) 

2.200 

Input  Cavity 


Drive  frequency  (GHz) 

5.5 

Operating  mode 

TM310 

Resistive  Quality  factor  (Q) 

2170 

Radius  (cm) 

1.701 

Length  (cm) 

0.225 

Magnetic  Field  Parameters 


Buncher  magnetic  field  (G) 

-567 

Output  magnetic  field  (G) 

613 

Cyclotron  frequency  (GHz) 

1.577 

Cusp  width  (cm) 

0.375 

Cavity  -  cusp  spacing  (cm) 

30 

Output  field  taper 

13% 

Output  cavity 


Output  frequency  (GHz) 

5.5 

Operating  mode 

TT-mode 

Diffractive  Q 

1240 

Resistive  Q 

13,430 

Length  (cm) 

15.000 

Radius  (cm) 

1.763 

Slot  depth  (cm) 

0.747 

Slot  width  (%) 

50 

Number  of  vanes 

6 
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0  7.231 

1  3.686 

2  5.055 

3  5.500 

Table  III.  Third  Harmonic  Amplifier  Performance 


Input  cavity  results 

Input  drive  power  (kW) 

1.16 

at  beam  (kV^m) 

11.3 

Drift  region  results 

Total  Av^  (%) 

9.07 

Total  Avx  (%) 

2.50 

Bunching  effectiveness  (%) 

37.79 

Energy  spread  A  E  (%) 

3.88 

Output  cavity  results 

Power  (kW) 

99.7 

Efficiency  (%) 

55.4 

Gain  (dB) 

19.3 

E^  at  beam  (kV/cm) 

21.3 
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Figure  Captions 


Figure  1 .  A  schematic  of  the  axially-modulated,  cusp-injected,  large-orbit  gyrotron. 

Figure  2.  Cross-sectional  view  of  the  vane  resonator  structure. 

Figure  3.  Input  cavity  field  radial  profiles:  axial  electric  field  (solid  line)  and  azimuthal  mag¬ 
netic  field  (dashed  line). 

Figure  4.  Output  cavity  field  radial  profiles:  azimuthal  electric  field  (solid  line)  and  axial  mag¬ 
netic  field  (dashed  line).  The  azimuthal  angle  is  fixed  at  the  center  of  a  vane  opening. 

Figure  5.  The  dependence  of  the  start-oscillation  current  on  magnetic  field  for  the  “7i”-mode. 

Figure  6.  Axial  evolution  of  the  radial  location  of  representative  particles  for  the  optimal  effi¬ 
ciency  parameters. 

Figure  7.  The  cross-sectional  distribution  of  representative  particles  for  the  optimal  efficiency 
parameters  at  (a)  the  entrance  of  the  output  cavity  and  (b)  the  exit  of  the  output  cav¬ 
ity. 

Figure  8.  The  dependence  of  output  efficiency  (solid  line)  and  total  axial  velocity  spread 
(dashed  line)  on  beam  thickness. 

Figure  9.  The  drive  curve  (solid  line)  for  the  zero-thickness  beam  when  the  output  cavity  Q  is 
~1 100.  The  dashed  line  indicates  the  bunching  efficiency. 

Figure  10.  The  dependence  of  output  efficiency  on  output  cavity  Q  for  the  zero-thickness  beam. 

Figure  1 1 .  The  dependence  of  output  efficiency  on  average  velocity  ratio  (solid  line).  The  re¬ 
quired  drift  tube  magnetic  field  is  indicated  by  the  dashed  line. 

Figure  12.  (a)  The  efficiency  as  a  function  of  harmonic  number  for  a  zero-thickness  beam  (solid 
line)  and  a  0.75  mm  thick  beam  (dashed  line).  The  output  cavity  length  is  fixed  at  15 
cm  and  the  magnetic  field  taper  is  fixed  at  13%,  but  the  drive  frequency  is  propor¬ 
tional  to  the  harmonic  number.  The  remaining  design  parameters  are  adjusted  to  be 
compatible  with  the  drive  frequency,  (b)  The  required  quality  factor  for  the  zero¬ 
thickness  beam  (solid  line)  and  the  finite  thickness  beam  (dashed  line). 


18 


Electron 

Gun 


Input  Drift  Pole 

Cavity  Region  Piece 


Output 

Cavity 


Figure  1.  A  schematic  of  the  axially-modulated,  cusp-injected,  large-orbit  gyrotron. 
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Figure  4.  Output  cavity  field  radial  profiles:  azimuthal  electric  field  (solid  line)  and  axial  mag¬ 
netic  field  (dashed  line).  The  azimuthal  angle  is  fixed  at  the  center  of  a  vane  opening. 
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Figure  6.  Axial  evolution  of  the  radial  location  of  representative  particles  for  the  optimal  effi¬ 
ciency  parameters. 
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Figure  7.  The  cross-sectional  distribution  of  representative  particles  for  the  optimal  efficiency 
parameters  at  (a)  the  entrance  of  the  output  cavity  and  (b)  the  exit  of  the  output  cav¬ 
ity. 
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Figure  10.  The  dependence  of  output  efficiency  on  output  cavity  Q  for  the  zero-thickness  beam. 
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Figure  12.  (a)  The  efficiency  as  a  function  of  harmonic  number  for  a  zero-thickness  beam  (solid 
line)  and  a  0.75  mm  thick  beam  (dashed  line).  The  output  cavity  length  is  fixed  at  15 
cm  and  the  magnetic  field  taper  is  fixed  at  13%,  but  the  drive  frequency  is  propor¬ 
tional  to  the  harmonic  number.  The  remaining  design  parameters  are  adjusted  to  be 
compatible  with  the  drive  frequency,  (b)  The  required  quality  factor  for  the  zero¬ 
thickness  beam  (solid  line)  and  the  finite  thickness  beam  (dashed  line). 
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On  the  design  of  electron  guns  for  cyclotron  resonance  masers 
utilizing  trochoidal  electron  beams  (trochotrons) 

W.  Lawson,'^  M.  Garven,*  and  G.  S.  Nusinovich 
Institute  for  Plasma  Research 
University  of  Maryland 
College  Park,  MD  20742 

Abstract 

The  properties  of  crossed-field  electron  guns  which  are  capable  of  providing 
the  required  trochoidal  beams  for  efficient  cyclotron  resonance  interaction  in  the 
trochotron  configuration  are  analyzed.  A  discussion  of  the  beam  parameters  that 
are  used  to  characterize  the  device  performance  is  included  in  a  general  descrip¬ 
tion  of  the  operation  of  these  guns.  Two  concrete  designs  are  presented  and  are 
used  to  demonstrate  the  performance  capabilities  of  these  devices  with  respect  to 
the  beam  self-fields,  the  externally  applied  electric  and  magnetic  fields,  and  the 
electrode  geometry.  It  is  shown  that  these  devices,  which  typically  generate  beams 
with  large  ratios  of  perpendicular  to  parallel  (drift)  velocity,  are  capable  of  pro¬ 
ducing  extremely  small  drift  velocity  spreads  over  a  wide  range  of  parameters. 


^Also  with  the  Electrical  Engineering  Department,  University  of  Maryland  at  College  Park. 
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L  Introduction 


It  is  well  known  that  a  non-relativistic  charged  particle,  in  the  presence  of  uniform,  perpen¬ 
dicular,  static,  external  electric  and  magnetic  fields,  will  execute  a  trochoidal  trajectory  which 
can  be  described  as  the  superposition  of  a  cyclotron  orbit  and  a  linear  drift.  This  drift  is  perpen¬ 
dicular  to  both  applied  fields  and  has  a  speed  equal  to  the  ratio  of  the  electric  field  magnitude  to 
the  magnetic  field  magnitude.  Cyclotron  resonance  masers  (CRMs),  based  on  the  interaction  of 
fast  electromagnetic  waves  with  these  trochoidal  beams  (trochotrons)  were  investigated  in  the 
early  days  of  CRM  development  [1,2]  and  have  recently  been  revisited  [3].  The  original  tro¬ 
chotrons  demonstrated  reasonably  good  efficiencies  in  traveling-wave  amplifier  configurations, 
~  but  were  later  abandoned  in  favor  of  CRM  devices  which  interacted  with  helical  beams 
(gyrotrons  [4-8]).  This  abandonment  occurred  because  the  initial  thrust  of  CRM  development 
was  for  millimeter  wave  operation  and  the  trochotron  has  an  inherent  frequency  limitation  that 
the  gyrotron  does  not.  Because  the  operating  frequency  in  a  CRM  device  is  tied  to  an  harmonic 
of  the  cyclotron  frequency,  the  applied  magnetic  field  must  increase  linearly  with  increasing  fre¬ 
quency.  And  since,  for  a  given  axial  drift  velocity,  the  electric  field  strength  in  a  trochotron  is 
tied  to  the  magnetic  field,  the  frequency  can  be  increased  only  until  the  breakdown  limit  is 
reached  in  the  device.  In  gyrotrons,  there  is  no  applied  dc  electric  field  in  the  interaction  region 
and  correspondingly  no  such  limitation. 

Still,  there  are  at  least  three  good  reasons  to  revisit  trochotron  operation.  First,  there  are  ap¬ 
plications  in  radar,  accelerators,  and  material  processing  that  require  centimeter-type  wave¬ 
lengths.  Second,  the  early  trochotron  experiments  were  designed  to  be  resonant  at  the  first 
harmonic  of  the  cyclotron  frequency.  Operation  of  CRM  devices  at  higher  harmonics  is  now 
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better  understood,  is  fairly  commonly  used  [9-11],  and  has  the  potential  to  increase  the  achiev¬ 
able  operating  frequency  of  trochotrons  significantly.  Finally,  trochotrons  can  potentially  be 
more  effective  in  devices  that  operate  fairly  far  from  the  cutoff  frequency  (e.g.  wide-band  ampli¬ 
fiers),  because  the  trochoidal  beam  has  considerably  less  axial  velocity  spread  than  the  helical 
beam.  So,  for  a  given  axial  wavenumber,  the  inhomogeneous  Doppler  broadening  of  the  cyclo¬ 
tron  resonance,  which  spoils  the  interaction  efficiency,  can  be  much  smaller  in  trochotrons  than 
in  gyrotrons.  This  is  true  because  the  axial  drift  velocity  is  tied  to  the  external  fields  and  the 
spread  arises  only  from  the  relative  space-charge  potential  depression  and  the  inhomogeneities  in 
the  external  fields.  In  helical  beams,  however,  the  axial  velocity  is  unconstrained  and  is  deter¬ 
mined  by  the  details  of  the  electron  gun  optics  and  the  space  charge  effects.  Typical  Magnetron 
Injection  Guns  (MIGs),  which  generate  beams  for  the  majority  of  gyrotrons,  theoretically  have 
perpendicular  velocity  spreads  of  a  few  percent.  [12-15]  In  addition,  experiments  with  high  cur¬ 
rent,  high  compression  MIGs  have  indicated  that  actual  perpendicular  spreads  may  be  signifi¬ 
cantly  higher.  [16,17]  Furthermore,  axial  velocity  spreads  in  helical  systems  are  typically  larger 
than  the  corresponding  perpendicular  spreads  by  the  square  of  the  average  perpendicular-to- 
parallel  velocity  ratio.  In  trochoidal  beams,  there  is  no  such  coupling  between  the  spreads  and 
extremely  large  ratios  of  perpendicular-to-parallel  velocities  (or  more  accurately  rotating-to-drift 
velocities)  can  be  achieved. 

At  the  University  of  Maryland,  we  have  begun  to  re-investigate  the  trochotron  configuration. 
In  an  earlier  paper  [3],  one  of  the  authors  analyzed  the  basic  geometry,  developed  some  scaling 
laws,  and  identified  some  important  restrictions  for  this  device.  It  was  pointed  out  that  the  suc- 
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cess  of  these  devices  is  contingent  on  the  ability  to  produce  a  high  quality  electron  beam  with  the 
required  parameters.  That  issue  is  addressed  in  this  paper. 

In  the  next  section  we  present  the  basic  configuration  of  a  planar  electron  gun  that  can  be 
used  to  generate  the  trochoidal  beam.  The  operation  of  the  device  is  described  and  the  method  of 
determining  the  beam  parameters  is  explained.  EGUN  [18]  simulations  of  a  simple  planar  gun 
geometry  are  performed  in  the  third  section  to  demonstrate  the  basic  properties  of  the  device.  A 
more  sophisticated  device  that  uses  an  injection  scheme  which  is  similar  to  an  earlier  Russian 
design  [1]  is  described  in  the  fourth  section.  Both  designs  are  geared  to  produce  the  necessary 
beam  parameters  for  a  moderate  power,  pulsed,  35  GHz,  second  harmonic  trochotron.  In  the  final 
section  of  this  paper  we  summarize  the  results  of  this  study. 

II .  The  basic  planar  electron  gun  geometry 

The  geometry  of  a  simple  planar  trochotron  electron  gun  is  shown  schematically  in  Figure  1. 
The  device  is  assumed  to  be  infinite  in  the  direction  perpendicular  to  the  figure  surface  so  that 
edge  effects  are  neglected.  The  cathode  is  centered  at  the  origin  and  the  distance  from  the  cath¬ 
ode  center  to  the  end  of  the  gun  is  /g.  The  emitter  length  is  4  and  the  cathode  loading  is  assumed 
to  be  uniform  with  a  density  of  /<.•  Thus,  the  total  beam  current  per  unit  length  is  4  =  4.  The 

gap  between  the  anode  and  cathode  at  the  emitter  strip  is  and  the  gap  between  the  plates  at  the 
beam  exit  is  d^.  The  magnetic  field  is  at  the  emitter  strip  and  at  the  end  of  the  gun.  The 
magnetic  compression  is  the  ratio  of  the  two  values: 

fm  =  BJB^.  (1) 

The  anode  -  cathode  potential  difference  is  V^.  Far  from  the  transition  region,  the  external 
electric  field  between  the  plates  is  uniform,  normal  to  the  electrode  surfaces,  and  equal  to  the 
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potential  difference  divided  by  the  gap  distance.  Thus,  E^=VJ  in  the  circuit  region  and  = 
Vq  !  th®  emitter  strip.  These  results  assume,  of  course,  that  the  electron  beam’s  self-fields 
can  be  neglected.  The  electric  compression  is  defined  by  the  ratio: 

f^  =  EJ  E^  =  d,/do.  (2) 

We  assume  that  electrons  are  emitted  from  the  cathode  with  virtually  no  kinetic  energy  (e.g. 
thermionic  emission).  Since  the  net  velocity  of  each  electron  is  the  sum  of  its  drift  velocity,  Vj, 
and  its  cyclotron  velocity,  v,.,  the  two  components  must  have  equal  magnitudes  but  opposite  di¬ 
rections.  Thus,  the  velocity  ratio  near  the  emitting  surface  is: 

Vc(0) /Vd(0)  =\,  (3) 

independent  of  the  gun  geometry.  The  velocity  ratio  at  the  end  of  the  gun  can  be  determined 
analytically  with  the  single  additional  assumption  that  the  magnetic  field  variation  is  adiabatic.  If 
we  define  a  to  be  this  ratio,  then 

^  =  /Vd(l^  =  [v,(0) /Vd(0)]  X [v,(ls)  /vc(0)] / [VdOs)  /^d(0)].  (4) 

The  first  term  in  (4)  is  just  equal  to  one  from  (3)  and  the  third  term  is  just 

VdOs)/Vd(0)  =  (E,/BJ/(E,/BJ  =/,//„  (5) 

firom  (1)  and  (2).  Under  adiabatic  conditions,  the  second  term  in  (4)  for  a  weakly  relativistic 
electron  beam  is  just  equal  to  the  square  root  of  the  magnetic  compression,  [14]  so 

^  ~  J  m  '  fe  •  (6) 

Another  parameter  of  primary  importance  is  the  average  energy  of  the  electrons  (eVb)  in  the 
gun  relative  to  the  applied  voltage  Vg.  A  non-relativistic  approximation  of  this  value  can  be 
quickly  determined  as  follows.  Neglecting  self-fields,  the  non-relativistic  particle  in  the  homoge- 
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neous  field  region  attains  its  average  energy  when  the  two  velocity  components  are  perpendicular 
to  each  other: 


2.  I 


From  (7),  the  relative  potential  can  be  readily  found: 

(8) 

JedlBl 

The  remaining  design  values  and  equations  come  from  the  consideration  of  the  microwave 
interaction  and  the  technical  limits  on  breakdown.  Bg  is  determined  by  the  resonance  condition 
and  Vq  and  Ip  are  derived  from  power  requirements.  The  spacing  d^  is  set  by  the  microwave 
mode,  given  the  realistic  constraints  on  Eq.  There  are  at  least  two  reasons  to  consider  a  lower 
value  of  electric  field  in  the  cathode  region  in  spite  of  the  fact  that  it  increases  the  required  mag¬ 
netic  compression  for  a  given  velocity  ratio.  The  first  is  that  the  restriction  on  the  maximum  al¬ 
lowable  electric  field  is  typically  more  stringent  near  the  emitter.  The  second  advantage  is  that 
the  relative  width  of  the  beam  compared  to  the  distance  it  travels  perpendicular  to  the  emitter  is 
smaller,  minimizing  the  effects  of  non-adiabatic  changes  in  the  electric  and  magnetic  field.  This 
in  turn  helps  to  minimize  the  spread  in  the  average  beam  energy.  A  potential  disadvantage  of  the 
lower  field  is  the  associated  decrease  in  the  space-charge  limited  current  density,  which  is 
given  for  the  infinite  planar  case  by  [19]: 


J  «  2.34  X  /  c/,' 


where  the  units  for  Jjc  V^,  and  d^  are  kicm,  V,  and  cm,  respectively. 


6 


In  the  following  two  sections,  we  will  use  these  equations  to  design  and  analyze  two  electron 
gun  configurations  which  are  suitable  for  a  35  GHz  tube.  The  two-dimensional  code  EGUN  is 
used  to  simulate  the  electron  guns,  so  edge  effects  are  neglected.  Convergence  is  checked  by 
varying  the  mesh  size  and  the  number  of  rays.  Magnetic  field  tapering  is  modeled  along  the  axis 
either  by  the  smooth  function: 


=  (B„  ^  \anh[(z  -  z„, )  /  ^ 


(10) 


or  by  the  piecewise  linear  function: 


°c  z  <  z„,  -  ^ 

\ 

'  2  2 

z>  z„.  +  ^ 

where  z„  is  the  midpoint  of  the  field  transition  and  ^  is  a  measure  of  the  transition  width.  These 
models  are  adequate  because  of  the  relative  insensitivity  of  the  beam  parameters  to  the  field 
profile.  To  model  inhomogeneities  in  the  applied  magnetic  field  of  a  real  magnetic  coil  system, 
we  use  (10)  and  stop  the  simulations  before  5,  attains  its  maximum  value. 

Characterization  of  the  average  velocity  ratio  (and  other  beam  parameters)  is  complicated  by 
space  charge,  inhomogeneities  in  the  applied  fields,  and  the  fact  that  the  relativistic  correction  to 
the  electron  motion  in  a  reference  frame  moving  at  the  drift  velocity  Vj  results  in  non-circular 
motion.  The  values  quoted  in  this  paper  are  found  by  averaging  the  simulated  quantities  at  vari¬ 
ous  locations  in  the  final  orbit  of  each  electron.  Uncertainties  in  these  quantities  are  always 
checked  and  are  typically  quite  small. 


z  -  z.. 


(11) 
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III.  An  in-line  gun  design 

The  electron  guns  characterized  in  this  paper  are  designed  to  meet  the  needs  of  a  second 
harmonic  trochotron  amplifier  with  the  parameters  given  in  Table  I.  The  specified  magnetic  field 
represents  a  12.6%  detuning  from  the  resonant  field  required  for  35  GHz  operation  at  the  second 
harmonic.  We  assume  a  thermionic  emitter  and  restrict  the  cathode  loading  accordingly;  the  im¬ 
pact  of  using  a  field-emitter  array  (FEA)  cathode  will  be  discussed  briefly  in  the  final  section. 
The  beam  power  density  is  selected  to  demonstrate  useful  power  levels  with  a  fairly  narrow  sys¬ 
tem.  The  peak  electric  field  indicated  is  for  the  circuit  region  and  is  below  what  has  been 
achieved  in  short  pulse  (microsecond)  systems.  We  restrict  the  peak  field  in  the  anode-cathode 
~  region  to  half  of  this  value. 

The  optimized  electrode  field  configuration  and  the  magnetic  field  profile  for  the  in-line  gun 
design  are  shown  in  Fig.  2.  The  simulated  beam  trajectory  is  also  indicated  for  the  nominal  pa¬ 
rameters  given  in  Table  II.  The  anode  cathode  gap  is  reduced  from  0.84  cm  in  the  emitter  region 
to  0.42  cm  over  a  distance  of  1.9  cm.  The  emitter  length  is  0.4  cm  and  the  surface  is  curved  to 
minimize  the  spread  in  the  beam  energy.  The  left-most  point  on  the  emitter  surface  is  parallel  to 
the  anode  and  the  surface  changes  by  10°  with  a  2.29  cm  radius  of  curvature.  The  current  density 
is  selected  to  achieve  the  design  beam  power  density  if  the  beam  is  centered  between  the  two 
plates.  The  cathode  is  shown  as  two  separate  conductors.  This  is  done  so  that  the  right  electrode 
could  in  principle  be  (slightly)  reversed  biased  to  minimize  the  production  of  secondary  elec¬ 
trons.  The  two  electrodes  were  kept  at  the  same  potential  in  all  of  the  simulations  that  we  per¬ 
formed. 
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As  evident  in  Fig.  2,  the  beam  is  well-formed  and  maintains  a  high  degree  of  phase  coher¬ 
ence  throughout  the  gun  region.  The  simulated  gun  parameters  at  the  nominal  design  point  are 
given  in  the  middle  column  of  Table  III.  The  average  velocity  ratio  of  5.0  is  obtained  with  a  drift 
velocity  spread  of  less  than  0.1%.  The  combination  of  these  two  values  is  unheard  of  in  MIG  de¬ 
signs  and  is  a  direct  consequence  of  the  drift  velocity  mechanism.  The  spread  in  pitch  angle, 
though  over  an  order  of  magnitude  higher,  is  still  considerably  smaller  than  what  is  readily  at¬ 
tainable  in  MIGs.  The  cathode  -  average  beam  position  spacing  is  about  56%  of  the  cathode  -  an¬ 
ode  spacing,  so  that  the  targeted  beam  power  is  modestly  exceeded.  The  beam  voltage  spread 
arises  in  part  from  self-fields,  but  mostly  from  the  2.3%  spread  in  the  electron  guiding  center 
radii.  The  simulated  space-charge  limiting  (SCL)  current  density  is  about  52  AJcm  and  is  higher 
than  the  theoretical  value  of  37  AJ crn  obtained  from  (9)  due  to  the  finite  emitter  width.  Clearly 
the  nominal  loading  of  10  AJcm  is  well  into  the  temperature-limited  regime.  The  peak  electric 
field  of  129  kV/cm  occurs  on  the  anode  at  the  end  of  the  gun.  The  electric  field  at  the  correspond¬ 
ing  point  on  the  cathode  is  about  18  kV/cm  lower  because  of  self-field  effects. 

The  dependence  of  average  velocity  ratio  on  the  applied  cathode  magnetic  field  (and  hence 
on  the  magnetic  compression)  is  indicated  in  Fig.  3a  by  the  solid  line.  The  theoretical  curve 
given  by  (6)  is  indicated  in  the  figure  by  the  dash-dotted  line.  Since  the  cathode  is  curved,  the 
average  cathode  -  anode  distance  is  used  in  the  equation.  Agreement  between  simulation  and  the¬ 
ory  is  seen  to  be  excellent  over  the  entire  range  of  fields.  Cathode  magnetic  fields  of  1.15  kG  and 
below  result  in  partial  beam  interception  on  the  anode.  Fields  above  1.65  kG  result  in  partial 
beam  interception  on  the  cathode.  Still,  the  velocity  ratio  can  be  adjusted  from  about  3.5  to  5.5  in 
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the  accessible  range  of  fields.  The  resultant  drift  velocity  spread  is  indicated  in  the  figure  by  the 
dashed  curve.  The  spread  remains  Avell  below  1%  until  the  field  exceeds  1.5  kG. 


From  (8),  it  is  evident  that  the  average  beam  energy  is  also  quite  dependent  on  the  cathode 
magnetic  field  (via  the  dependence  of  a  on  f^,).  The  simulated  and  analytic  dependencies  of  beam 
voltage  on  cathode  magnetic  field  are  indicated  in  Fig.  3b  by  the  solid  and  dash-dotted  lines,  re¬ 
spectively.  The  two  curves  indicate  a  similar  drop-off  in  voltage  with  decreased  magnetic  com¬ 
pression,  but  the  theoretical  curve  is  consistently  about  2  -  3  kV  higher.  At  least  half  of  this 
discrepancy  can  be  accounted  for  by  the  potential  depression  of  the  beam.  The  spread  in  beam 
voltage  is  shown  in  Fig.  3b  by  the  dashed  line.  As  with  the  spread  in  a,  there  is  a  marked  in¬ 
crease  after  the  field  exceeds  1 .5  kG.  Comparing  the  results  in  the  figure,  it  is  clear  that  the  pre¬ 
eminent  effect  of  cathode  magnetic  field  variations  is  on  the  beam  voltage  average  and  spread. 

The  dependence  of  average  velocity  ratio  on  the  voltage  applied  across  the  two  electrodes  is 
indicated  by  the  solid  line  in  Fig.  4.  Even  though  there  is  an  increase  in  the  electric  field  every¬ 
where  in  the  circuit,  the  electric  compression  remains  unchanged  and  (6)  predicts  no  variation  in 
a.  The  slight  increase  in  a  with  increasing  voltage  is  mainly  due  to  the  relative  decrease  in  the 
potential  depression.  Partial  beam  interception  on  the  anode  occurs  with  an  applied  voltage  of  60 
kV.  It  is  interesting  to  note  that  the  ratio  of  average  beam  voltage  to  applied  voltage  increases 
from  37%  at  40  kV  to  68%  at  60  kV.  The  drift  velocity  spread  is  given  in  the  figure  by  the 
dashed  line.  The  spread  is  well  below  0.5%  for  all  beam  voltages  above  45  kV.  The  decrease  in 
beam  quality  below  that  point  is  due  in  large  part  to  the  drop  off  in  beam  voltage  and  the  resul¬ 
tant  increase  in  relative  space-charge  forces. 
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Figure  5  demonstrates  the  dependence  of  average  velocity  ratio  (solid  line)  and  drift  velocity 
spread  (dashed  line)  on  cathode  loading.  Again,  the  simple  theory  predicts  no  variation  in  veloc¬ 
ity  ratio  and  the  weak  dependence  exhibited  in  the  simulations  can  be  directly  attributed  to  space- 
charge  depression.  This  small  change  in  a  could  readily  be  compensated  for  by  a  small  change  in 
the  cathode  magnetic  field,  if  desired.  A  high  quality  beam  is  achieved  from  current  densities  of 
2  A/cm  (J5  —  5  A/cm  )  up  to  5  A/cm  (J^  =  12.5  A/cm^).  The  increase  in  velocity  spread  with  in¬ 
creasing  current  is  due,  of  course,  to  space-charge  effects.  The  smaller  increase  in  spread  that  oc¬ 
curs  as  the  current  density  goes  to  zero  results  form  the  inhomogeneities  in  the  magnetic  field 
and  the  large  spread  in  guiding  center  transverse  locations.  This  spread  can  be  reduced  by  opti¬ 
mizing  the  electrode  surfaces  at  lower  current  densities,  but  that  effort  would  result  in  increased 
spread  at  the  operating  current  density. 

IV.  An  injection  gun  design 

An  injection  gun  design,  based  on  the  circuit  parameters  in  Table  I  but  utilizing  a  geometry 
similar  to  the  earlier  Russian  design  in  [1],  is  shown  in  Fig.  6.  Such  a  configuration  can  be  used 
to  isolate  the  microwave  path  from  the  beam  source  in  an  amplifier  configuration.  In  addition  to 
the  electrode  configuration,  the  magnetic  field  profile  and  the  simulated  ray  trajectories  are  also 
indicated  in  the  figure.  The  piecewise  linear  magnetic  field  profile  in  (11)  is  used  for  this  design 
because  of  the  fairly  large  transverse  range  of  the  beam.  The  shape  of  the  electrodes  in  the  emit¬ 
ter  region  is  nearly  identical  to  that  of  the  in-line  gun.  In  fact,  the  electrode  parameters  are  identi¬ 
cal  to  those  listed  in  Table  II.  Exactly  the  same  emitter  shape  and  nominal  cathode  loading  are 
used.  The  anode  has  the  same  14°  slope  as  the  in-line  design,  beginning  at  the  same  relative  axial 
position.  However,  the  taper  ends  after  about  5  mm  in  order  to  begin  the  transition  to  the  upper 
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plate.  To  promote  an  adiabatic  transition,  the  radius  of  the  curvature  of  the  lower  plate  in  the 
transition  region  is  selected  to  be  about  10  times  larger  than  the  final  Larmor  radius.  The  circuit 
magnetic  field  (Bq)  is  identical  in  the  two  cases,  but  the  complex  injection  geometry  requires  a 
cathode  field  of  Bg=1.355  kG  to  achieve  the  desired  velocity  ratio  of  5.0.  The  transition  width  is 
4=1 .9  cm  and  the  transition  center  is  z^=2A  cm. 

The  simulated  results  for  the  configuration  are  given  in  the  final  column  of  Table  III.  In  spite 
of  the  more  complex  geometry,  the  spread  in  the  drift  velocity  is  only  about  0.1%  and  is  still  well 
over  an  order  of  magnitude  smaller  than  the  best  values  found  in  MIGs.  The  average  beam  power 
is  only  slightly  lower  in  the  injection  design  because  of  a  small  difference  in  the  average  guiding 
center  location.  However,  the  guiding  center  spread  is  smaller  for  this  case,  and  this  results  in 
improved  values  of  beam  voltage  and  velocity  ratio  spread.  Comparing  the  beam  trajectories 
shown  in  Figs.  2  and  5,  it  is  apparent  that  the  injection  gun  beam  is  considerably  more  mixed 
than  the  in-line  beam.  The  mixing  begins  in  the  injection  region  but  doesn’t  appear  to  have  any 
adverse  effects  on  the  beam  quality.  The  peak  electric  field  occurs  on  the  lower  electrode  at  about 
3.5  cm  (in  the  middle  of  the  small  radiused  transition).  The  peak  field  on  the  anode  is  135  kV/cm 
and  occurs  on  the  radiused  end.  The  latter  value  may  be  adequate  and  the  higher  value  could  be 
reduced  in  an  actual  design  simply  by  increasing  the  transition  radius. 

V.  Discussion 

This  study  has  demonstrated  the  feasibility  of  electron  guns  for  trochotrons  in  a  parameter 
regime  of  interest  for  current  applications.  The  beam  voltage,  peak  electric  field,  and  cathode 
loading  are  all  consistent  with  what  is  readily  achievable  in  modern  tubes.  Ratios  of  transverse 
velocity  to  drift  velocity  of  5  (or  higher)  are  easily  obtainable  with  drift  velocity  spreads  well 
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below  1  Electron  beams  with  such  low  drift  velocity  spreads  are  definitely  suitable  for  reso¬ 
nant  operation  far  from  cutoff.  So,  the  main  issue  in  the  development  of  wideband  traveling- 
wave  trochotrons  could  be  providing  proper  coupling  between  the  electrons  and  a  wave  propagat¬ 
ing  with  a  large  group  velocity  (rather  than  the  deleterious  effect  of  inhomogeneous  Doppler 
broadening  inherent  in  gyro-TWTs). 

The  simple  analytic  model  of  the  gun’s  beam  dynamics  results  in  good  estimates  of  the  sys¬ 
tem  parameters  as  long  as  the  electric  and  magnetic  fields  are  varied  adiabatically  and  the  relative 
space-charge  effects  are  sufficiently  small.  Beam  quality  is  seen  to  remain  high  over  a  wide 
range  of  parameters.  Frequently,  the  parameter  most  sensitive  to  system  variations  is  the  relative 
location  of  the  beam  between  the  plates  and  the  resulting  average  beam  energy. 

Though  not  presented  here,  we  have  also  investigated  the  performance  of  trochotron  guns 
with  narrower  emitters  and  higher  current  densities.  Such  emitters  are  consistent  with  the  FEA 
cathodes  that  are  under  development  at  many  locations.  [20-22]  These  systems  typically  require 
no  curvature  and  give  higher  beam  quality  over  a  greater  range  in  current  densities. 
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Table  I.  The  microwave  circuit  requirements. 


Parameter 

Value 

Output  frequency  (GHz) 

35 

Applied  voltage  (kV) 

50 

Beam  power  density  (kW/cm) 

100 

Circuit  magnetic  field  (kG) 

6.0 

Maximum  electric  field  (kV/cm) 

150 

Nominal  cathode  loading  (A/cm^) 

10 

Table  II.  The  in-line  gun  parameters. 


Parameter 

Value 

Circuit  plate  gap  do  (cm) 

0.42 

Cathode-anode  plate  gap  d,,  (cm) 

0.84 

Nominal  Current  density  (A/cm) 

4.0 

Cathode  magnetic  field  B^  (kG) 

1.26 

Field  transition  location  (cm) 

1.45 

Field  transition  length  ^  (cm) 

0.6 

Cathode  curvature  (degrees) 

10 

Anode  taper  angle  (degrees) 

14 

14 


Table  III.  The  simulated  properties  of  the  two  gun  designs. 


Parameter 

In-Line 

Injection 

Average  velocity  ratio  a 

5.00 

5.00 

Velocity  ratio  spread  (%) 

1.70 

0.74 

Normalized  average  drift  velocity 

0.0607 

0.0605 

Drift  velocity  spread  (%) 

0.06 

0.11 

Average  beam  power  (kW/cm) 

106.4 

104.7 

Beam  voltage  spread  (%) 

3.55 

1.50 

SCL  current  density  (A/cm  ) 

52 

52 

Peak  electric  field  (kV/cm) 

129 

153 

15 
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Figure  Captions 

Figure  1 .  Layout  of  a  simple  trochoidal  gun  configuration. 

Figure  2.  Schematic  of  the  in-line  gun  electrode  configuration.  The  simulated  ray  trajecto¬ 
ries  at  the  nominal  parameters  of  Table  II  are  shown  in  the  figure  along  with  the 
magnetic  field  profile. 

Figure  3.  (a)  The  dependence  of  simulated  average  velocity  ratio  (solid  line)  and  drift  ve¬ 
locity  spread  (dashed  line)  on  cathode  magnetic  field.  The  theoretical  curve  for  the 
average  velocity  ratio  based  on  (6)  is  indicated  by  the  dash-dotted  line,  (b)  The 
dependence  of  simulated  average  beam  voltage  (solid  line)  and  voltage  spread 
(dashed  line)  on  cathode  magnetic  field.  The  theoretical  curve  for  the  average 
voltage  based  on  (8)  is  indicated  by  the  dash-dotted  line. 

Figure  4.  The  dependence  of  simulated  average  velocity  ratio  (solid  line)  and  drift  velocity 
spread  (dashed  line)  on  applied  voltage. 

Figure  5.  The  simulated  average  velocity  ratio  (solid  line)  and  drift  velocity  spread  (dashed 
line)  as  a  function  of  beam  current  density. 

Figure  6.  Schematic  of  the  injection  gun  electrode  configuration.  The  simulated  ray  trajec¬ 
tories  at  the  nominal  parameters  are  shown  in  the  figure  along  with  the  magnetic 
field  profile. 
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Figure  1 .  Layout  of  a  simple  trochoidal  gun  configuration. 
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Figure  2.  Schematic  of  the  in-line  gun  electrode  configuration.  The  simulated  ray  trajecto¬ 
ries  at  the  nominal  parameters  of  Table  II  are  shown  in  the  figure  along  with  the 
magnetic  field  profile. 
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Figure  3.  (a)  The  dependence  of  simulated  average  velocity  ratio  (solid  line)  and  drift  ve¬ 
locity  spread  (dashed  line)  on  cathode  magnetic  field.  The  theoretical  curve  for  the 
average  velocity  ratio  based  on  (6)  is  indicated  by  the  dash-dotted  line,  (b)  The 
dependence  of  simulated  average  beam  voltage  (solid  line)  and  voltage  spread 
(dashed  line)  on  cathode  magnetic  field.  The  theoretical  curve  for  the  average 
voltage  based  on  (8)  is  indicated  by  the  dash-dotted  line. 
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Figure  4.  The  dependence  of  simulated  average  velocity  ratio  (solid  line)  and  drift  velocity 
spread  (dashed  line)  on  applied  voltage. 
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Figure  5.  The  simulated  average  velocity  ratio  (solid  line)  and  drift  velocity  spread  (dashed 
line)  as  a  function  of  beam  current  density. 
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Figure  6.  Schematic  of  the  injection  gun  electrode  configuration.  The  simulated  ray  trajec¬ 
tories  at  the  nominal  parameters  are  shown  in  the  figure  along  with  the  magnetic 
field  profile. 
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